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PREFACE 

This  report  is  the  second  of  a  series  describing  experimental  work  presently 
under  way  at  Philco-Ford  Corporation  on  the  character  of  compressible  free 
turbulent  flows,  and  especially  on  the  characteristics  of  hypersonic  wakes. 
It  deals  with  the  measurement  and  the  spatial  and  spectral  resolution  of 
fluctuations  in  an  axi-symmetric  compressible  wake  and  discusses  the  sim¬ 
ilarity  rules  by  which  these  results  can  be  extended  to  higher  speeds. 

This  work  has  been  preceded  by  Philco-Ford  Report  No.  U-3978,  Mean-Flow 
Measurements  in  an  Axi-sy  letric  Compressible  Turbulent  Wake,  1  March  1967, 
in  which  measurements  of  the  mean  properties  of  the  wake  dealt  with  herein 
were  presented.  A  t'i  ird  report,  "Turbulent  Front  Structure  of  an  Axi 
symmetric  Compressible  Wake,"  describing  the  intermittency  measurements, 
is  being  prepared. 

The  author  acknowledges  with  gratitude,  the  financial  support  and  encourage¬ 
ment  of  the  Advanced  Research  Projects  Agency,  the  Space  and  Missile  Systems 
Organization,  USAF,  and  the  Philco-Ford  Corporation;  the  Analysis  Section  of 
the  latter  was  especially  helpful  with  the  machine  programming  and 
computations.  Special  thanks  are  also  due  Dr.  John  Laufer,  Consultant, 
for  many  stimulating  discussions,  and  Mr.  Lee  Von  Seggern,  whose  expert 
assistance  with  the  equipment  and  instrumentation  proved  of  great  value  In 
this  work. 


ABSTRACT 


The  turbulence  characteristics  of  an  axi-symmetric  we.ce  shed  by  a  very 
slender  body  at  Mach  3  have  been  measured  with  the  hot-wire  anemometer. 

The  region  covered  beg<vn  at  transition,  located  immediately  downstream  of 
the  model,  and  ended  about  100  virtual  wake  diameters  downstream.  By  using 
computer-aided  corrections  to  the  anemometer  frequency  response,  maps  of 
the  temperature  and  axial  velocity  fluctuations,  their  cross-correlations, 
spectral  density,  and  auto  correlation  macroscales  were  obtained  as  a 
function  of  the  axial  and  radial  coordinates.  Beyond  40  wake  diameters  from 
the  virtual  origin,  dynamic  equilibration  forces  these  distributions  into 
seemingly  self-preserving  forms.  These  forms  are  identical  with  their 
imcompressible  counterparts  except  for  a  coordinate  transformation  already 
suggested  by  the  dynamic  equilibrium  hypothesis  and  the  Reynolds  analogy 
typical  of  adiabatic  flows.  The  density  macroscale  is  numerically  smaller 
than  that  of  the  velocity,  although  they  both  approximate  the  transverse 
wake  scale.  Away  from  the  axis,  both  the  fluctuation  magnitudes  and  their 
spectral  densities  are  distorted  by  the  pseudoturbulence  generated  by  the 
intermittent  boundary,  which  was  found  to  be  a  weakly  periodic  structure 
with  a  wavelength  on  the  order  of  the  wake  diametpr. 
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Transfer  function  of  voltmeter  detector 
Dimensionless  hot-wire  overheat  (Equation  (21)) 
Transfer  function  of  amplifier  (Equation  (6)) 

Heat  capacity  of  hot-wire 

Function  of  hot-wire  finite  circuit  factor 

Virtual  turbulent  wake  diameter  (squared) 
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Ideal  mean-square  output  of  system  at  frequency  f 

Vacuum-tube  (mean-square)  voltmeter  output  of 
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Error  ratio  (Equation  (16)) 

Function  of  hot-wire  overheat 
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Hot-wire  length 
Transverse  scale  of  wake 
Local  mass-flux  (product  pu) 

Mach  number 

Wire  time  constant 

Hot-wire  constant  at  current  I 

Time  constant  of  amplifier 
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Hot-wire  Nusselt  number  based  on  stagnation 
conductivity 

Overall  transfer  function 

pressure 

Pitot  pressure 

Density  defect  (Equation  (86)) 

Hot-wire  resistance 

Mass-f lux-total  temperature  cross-correlation 
coefficient 

Velocity -temperature  cross-correlation  coefficient 
Hot-wire  re.  istance  at  current  I 
Recovery  resistance  of  wire' 

Heated  resistance  of  wire 
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Wire  Reynolds  number  (ba: ed  on  stagnation  viscosity) 
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Nondimensional  total  temperature  fluctuation 
Temperature 

Total  (local)  temperature 
Nondimensional  temperature  (Equation  (88)) 

Axial  velocity 

Nondimensional  axial  velocity  (Equation  (87)) 

Local  velocity  scale  (Equation  (57)) 

Wave  analyzer  output  (volts  dc) 

Same  as  V 

No’se  output  of  wave  analyzer  (volts  dc) 

Velocity  defect  (Equation  (64)) 


W  Compensator  transfer  function 

WEB  Designation  of  present  experiment 

WEB-Q  Q  (  ■  I,  II . )  designation  of  computer  program 

X  Axial  coordinate 

X  Nondimensional  axial  coordinate  (Equation  (81)) 

XQ  Axial  position  of  virtual  origin 

X-Station  Axial  position  of  measurements 

Radial  coordinate 

Function  of  Mach  number  (Equation  (37)) 
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Function  of  the  Mach  number  (Equation  (38)) 

Ratio  of  specific  heat 

Rms  value  of  the  local  variable  Q 

Frequency  interval 

Finite  circuit  factor 

Hot-wire  recovery  factor 

Nondimensional  radial  coordinate  (Equation  (60)) 

Temperature  defect  (Equation  (85)) 

Autocorrelation  macroscale  of  the  density 
fluctuations 

Autocorrelation  macroscale  of  the  velocity 
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Local  density  of  flow 

Density  scale 

Nondimensional  density  (Equation  (89)) 
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Same  as  above  but  within  a  passband  Af 
(i.e.,  density  fluctuation  density) 
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T  Rms  nondimensional  frequency-integrated  velocity 
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T(0  Same  as  above  but  within  a  passband  Af 

(i.e.,  velocity  fluctuation  density) 

Function  of  the  wire  overheat 
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SECTION  1 

INTRODUCTION 


Renewed  interest  in  the  study  of  turbulence  in  compressible  wakes  has  been 
generated  by  the  need  to  predict  the  radar-scattering  characteristics  of 
re-entry  flow  fields.  Formidable  difficulties  in  fluid  mechanics  and 
chemistry  have  directed  current  research  toward  two  distinct  goals.  One 
is  the  study  of  the  wake  turbulent  structure  itself;  the  other  is  a  search 
for  the  link  between  the  gasdynamical  and  electronic  turbulence  properties 
of  a  hydrodynamics lly  turbulent  plasma.  This  report  addresses  itself  to 
the  first  problem. 

Because  high-speed  wakes  involve  fluctuations  in  the  pressure,  density, 
and  temperature  in  addition  to  the  velocity,  unknowns  of  interest  should 
now  include  the  density  or  temperature  fluctuation  magnitude  and  its  radial 
and  axial  distributions,  as  well  as  the  distribution  of  spectra,  autocor¬ 
relation  functions,  and  scales  of  the  same  properties.  To  these,  one  should 
add  the  statistical  distribution  of  fluctuation  magnitudes  and,  of  course, 
the  intermittency  characteristics  of  the  turbulent  bcttndary.  Further,  there 
#  are  two  idiosyncracies  of  the  compressible  wake  which  bear  heavily  on  the 
design  of  an  experiment.  First,  the  inviscid  gradients  generated  by  vehicle 
shock  waves  force  the  wake  to  grow  into  a  nonuniform  stream;  this  puts  a 
premium  on  the  search  for  the  self-preserving  characteristics  of  the  turbulent 


wake  so  that  the  turbulent  structure  can  be  predicted  regardless  of  inviscid 
gradients.  Second,  in  contrast  to  low  speeds,  transition  in  the  supersonic 
wake  can  occur  at  some  distance  behind  the  vehicle;  the  question  arises  as 
to  the  initial  values  of  the  turbulence  properties  in  the  transition  zone. 

Very  little  has  been  done  to  answer  these  questions.  Most  of  the  available 
information  comes  from  two  sources;  low-speed  experiments  (References  1 
through  5)  and  very  recent  attempts  to  measure  these  properties  in  Lallistic 
ranges  (References  6  and  7).  The  attractiveness  of  producing  wakes  with 
hypervelocity  pellets  in  the  latter  is  severely  offset  by  two  factors; 
first,  the  lack  of  pellet  path  control  and  probe  position;  and,  second, 
difficulties  with  Taylor's  hypothesis  (Reference  8).  For  example,  it  is 
impossible  to  perform  repeated  shots  with  slender  bodies  at  zero  angle  of 
attack  and  to  control,  a  priori,  the  path  of  the  probe  in  the  wake;  the 
radial  position  of  the  measurement fpoint  is  therefore  always  in  doubt,  and 
measurements  near  the  wake  axis  obviously  cannot  be  made  for  many  tens 
(perhaps  hundreds)  of  base  diameters.  However,  at  these  far  distances, 
the  mean  flow  velocity  around  the  probe  might  be  »o  low  that  the  turbulence 
structure  might  change'  significantly;  the  so-called  Taylor's  hypothesis, 
which  is  crucial  to  the  interpretation  of  data,  is  thus  invalid. 

We  conclude  that  ballistic  range  experiments,  although  of  qualitative  value 
at  high  Mach  numbers,  cannot  match  the  precision  and  accuracy  of  wind  tunnel 
work  in  this  regard. 

In  the  meantime,  attempts  have  been  made  to  approach  the  problem  theoretically; 
several  are  summarized  in  Reference  9.  Recently,  Morkovin  (Reference  10) 
correlated  low-  and  high-speed  experimental  data  within  the  scope  of  the  so- 
called  dynamic  equilibrium  hypothesis  and  his  suggestions  were  expanded  by 
Laufer  (Reference  11)  specifically  to  the  axi-symmetric  hypersonic  wake. 

The  power  of  this  latter  approach  lies  in  the  dependence  on  concepts  verified 
experimentally  (Reference  1)  and  on  its  potential  of  providing  useful  pre¬ 
dictions  for  the  hypersonic  wake  by  unifying  incompressible  and  compressible 
behavior  into  relatively  simple  rules.  By  virtue  of  both  its  simplicity  and 

promise,  this  theory  will  play  a  central  role  in  interpreting  the  present 
experimental  data. 
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No  wind  tunnel  turbulence  measurements  of  compressible  wake  flows  have  been 
reported  to  date.  Such  experiments  are  held  back  mainly  by  shortcomings  of 
hot-wire  anemometry  and,  in  the  case  of  axi-symmetric  wakes,  by  the  diffi¬ 
culty  of  suspending  axi-symmetric  models  at  high  speeds.  Of  these  two 
problems,  the  former  is  the  most  serious,  and  was  overcome  in  the  present 
experiment  by  a  complex  and  admittedly  laborious  process. 

The  present  work  deals  with  the  turbulent  field  alone  and  has  been  preceded 
by  a  detailed  study  of  the  mean  (time-averaged)  flow  field  in  the  axi- 
symmetric  wake.  These  mean-flow  measurements,  described  in  Reference  12, 
are  necessary  in  order  to  compute  the  sensitivity  of  the  hot  wire  to  the 
different  modes  of  turbulence  and  to  relate  the  turbulence  properties  to 
the  local  and  integral  properties  of  the  flow.  The  reader  should  consult 
Reference  12  for  details  regarding  the  mean-flow  properties  and  the  instru¬ 
mentation  and  techniques  by  which  these,  properties  were  measured. 


SECTION  2 


APPARATUS  AND  TECHNIQUE 


2.1  WIND  TUNNEL  AND  MODEL 

The  experiment  was  performed  at  Mach  3  In  a  continuous  air  wind  tunnel. 

The  wake  was  produced  by  a  circular  rod  supported  by  a  strut  upstream  of 
the  nozzle  throat  and  cantilevered  into  the  test  section  through  the 
throat.  The  Reynolds  number  in  the  test  section  was  50,000  per  cm, 
generated  by  a  tunnel  total  pressure  and  temperature  of  508  mm  Hg  abs 
ard  27°C,  respectively.  Further  details  on  the  tunnel  itself  appear  in 
Reference  13. 

2.2  GENERAL  FEATURES  OF  THE  FLOW  FIELD 

Figure  1  shows  a  schematic  view  of  the  model  and  wake  flow.  The  wake  is 
formed  by  the  boundary  layer  shed  off  the  body;  transition  to  turbulence 
could  be  located  will  along  the  axi-symmetric  wake  by  changing  the 
tunnel  total  (supply)  pressure.  At  the  chosen  pressure  of  508  ram  Jig, 
transition  Ives  between  about  10  and  15  base  diameters  downstream  of  the 
base,  allowing  for  the  maximum  turbulent  wake  length  accessible  to  probes; 
at  the  same  time,  the  wake  is  relatively  thinner  than  would  be  obtained  if 
transition  was  moved  to  the  model  surface.  Furthermore,  since  laminar 
reattachment  (at  the  wake  neck)  incurs  a  lower  pressure  rise  than  turbulent 
reattachment,  inviscid  gradients  around  the  wake  were  extremely  mild;  the 


FIGURE  1. 


TOP:  PHOTOGRAPH  OF  TEST  SECTION  SHOWING  AXISYMMETRH 

PR0BE  AT  LEFT  AND  PROBE  ACTUATOR  (TOP  LEFT). 
BOTTOM:  DETAILS  OF  THE  AXISYMMETRIC  MODEL.  ' 
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wake  shock  disappeared  completely  at  this  pressure  and  the  Mach  number 
external  to  the  turbulent  region  was  2.96  ±0.02  for  the  entire  wake  length. 

2.3  HOT-WIRE  ANEMOMETER 

The  turbulence  properties  were  measured  with  a  single  0.00005-inch-diameter, 
0.01-inch-long  pt  107.  Rh  wire;  the  hot-wire  probe  is  pictured  in  Figure  2. 
Prior  to  use,  each  wire  underwent  a  temperature-resistance  calibration  in 
a  controlled  oven  and  a  calibration  of  its  heat-transfer  characteristics 
in  a  known  supersonic  flow  external  to  the  wake.  The  details  of  these 
preparatory  steps  are  given  in  Reference  12,  together  with  numerical 
results  of  the  calibrations. 

The  single  most  important  precaution  in  mounting  the  hot  wire  for  use  was 
its  strain  gage  characteristic.  The  problem,  arising  from  structural 
oscillations  of  the  wire  in  the  flow  and  superposing  on  the  output  spectrum 
numerous  misleading  energy  peaks,  has  been  thoroughly  described  by  Morkovin 
(Reference  14).  By  using  only  those  wires  with  significant  slack  (e.g.,  a 
quarter-circle),  it  was  possible  to  obtain  completely  undistorted  spectra. 
In  practice,  each  single  wire  was  routinely  subjected  to  a  test  whereby  its 
output  spectrum  in  the  turbulent  wake  was  obtained  prior  to  use.  A  wave 
analyzer  with  a  mechanically  driven  frequency  sweep  was  especially  helpful 
in  these  tests.  In  practice,  about  80  percent  of  the  wires  mounted  were 
rejected  because  the  spectrum  so  obtained  showed  discernible  strain  gage 
problems ,  of  the  type  shown  in  Figure  3. 

Spatial  resolution  of  the  wire  in  the  transverse  direction  was  controlled, 
of  course,  by  its  length,  which  was  0.01  inch.  By  contrast,  the  wake 
diameter  was  about  0.30  inch  on  the  average  so  that  the  probe  resolution 
was  of  order  1/30.  It  should  be  also  noted  that  the  wire  was  positioned 
at  all  times  along  a  direction  normal  to  the  distance  between  it  and  the 
wake  axis.  In  addition  to  the  spatial  resolution  necessary  for  these 
measurements,  the  temporal  resolution  was  of  great  importance.  This  was 
controlled  by  the  wire  time-constant  limit  at  zero  current  (on  the  order 
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FIGURE  3.  STRAIN -GAGE  CHARACTERISTICS  OF  HOT-WIRES  PLOTTED  AS  INTEGRATED 
WIRE  OUTPUT  VS.  FREQUENCY.  TOP  LEFT:  SEVERE  TOP  RIGHT: 

TYPICAL  BOTTOM  LEFT:  ACCEPTABLE  BOTTOM  RIGHT:  IDEAL  (CURVES 
RL PRODUCED  DIRECTLY  FROM  i-Y  PLOTTER  SPECTRA) 
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of  0.1  millisecond)  and  the  f loor-to-ceiling  ratio  of  the  compensating 
amplifier,  which  was  500  in  the  present  instance.  The  resulting  frequency 
response  of  800  kilocycles  per  second  represented,  when  divided  by  the 
free-stream  velocity,  a  scale  of  about  0.2  the  body  diameter.  This  short¬ 
coming  was  significant  but  not  prohibitive  and,  as  will  be  seen  below,  was 
overcome  by  detailed  spectral  measurements  at  the  higher  frequencies.  The 
frequency  response  of  the  anemometer  (its  transfer  function)  is  discussed 
in  Appendix  A.  An  idea  of  the  response  speed  of  the  instrument  under 
typical  operating  conditions  can  be  obtained  from  Figure  4. 


2.4  ELECTRONIC  INSTRUMENTATION 

A  block  diagram  of  the  instrumentation  necessary  to  record  the  turbulence 

data  is  shown  in  Figure  5;  the  same  equipment  was  earlier  used  for  the 

* 

mean  measurements. 


The  wire  was  powered  by  mercury  cells  contained  in  the  Transmetrics,  Inc. 
Model  6401-1  Current  Control  Panel,  and  signal  amplification  and  compensa¬ 
tion  was  achieved  by  a  Transmetrics  Model  6401-5  Constant  Current  Amplifier 
and  a  Model  6401-7  Square  Wave  Generator.  A  Hewlett-Packard  Model  3440A 
Digital  Voltmeter  and  a  switching  circuit  selectively  displayed  the  wire 
current  and  voltage,  alleviating  the  need  of  bridge -balancing  during  the 
probe  traverse.  Measurement  and  display  instruments  included  a  Ballantine 
Model  320A  true  RMS  Voltmeter,  a  Singer-Metrics  Model  SB-76Z  1  to 
300-kilocycle-per-second  Panoramic  Wave  Analyzer,  a  Hewlett-Packard 
Model  310A  1.5-megacycle  Wave  Analyzer  driven  by  a  297A  Sweep  Drive,  a 
Tektronix  Model  551  Dual-Beam  Oscilloscope,  and  a  Hewlett-Packard 
Model  7000A  X-Y  Recorder.  An  intermit tency  meter,  designed  by 
Dr.  H.  Shapiro  and  built  at  Philco-Ford  by  Mr.  L.  Von  Seggern,  was  used 
for  intermittency  and  other  statistical  measurements.  A  Tektronix  Model  549 
Storage  Oscilloscope  with  dual-trace  capability  (effected  by  a  Tektronix 
Model  1A1  Chopper)  was  also  used  for  statistical  measurements. 
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FIGURE  4.  TYPICAL  OSCILLOGRAM  OF  INTERMITTENT  WIRE  OUTPUT  NEAR  WAKE  EDGE. 

TOP  TRACE  IS  VOLTAGE  OUTPUT  SHOWING  A  15-MICROSECOND  "EVENT' 
DURING  WHICH  THE  WIRE  IS  IMMERSED  IN  TURBULENCE.  BOTTOM  TRACE 
IS  RECTIFIED  ANALOG  OF  THE  TOP  TRACE,  USED  FOR  MEASURING  THE 
INTERI  ITTENCY  FACTOR.  NOTE  FAST  HOT-WIRE  RESPONSE:  DOTS  ARE 
SPACED  ONE  MICROSECOND  APART. 
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FIGURE  5.  ELECTRONIC  CIRCUIT  USED  FOR  TURBULENCE  MEASUREMENTS 


2 . 5  EXPERIMENTAL  PROCEDURE 


The  measurement  of  the  rms  fluctuations  of  the  flow  variables  and  of  their 
spatial  ccrrelations  in  the  wake  lies  within  the  capability,  of  the  constant 
current  anemometric  technique  developed  for  compressible  flows  by  Kovasznay 
Morkovin,  and  their  co-workers  (References  14  and  15).  The  latter  tech¬ 
nique,  bolstered  by  the  computer-aided  process  of  extending  the  wire 
frequency  response,  was  used  throughout  this  work.  The  measurements  were 
performed  by  obtaining  radial  profiles  of  hot-wire  output  at  14  axial 
positions  along  the  wake  (designated  X-Station  0  through  13)  beginning  at 
17  model  base  diameters  downstream  of  the  base  and  spaced  0.500  inch  apart. 
The  several  weeks  needed  to  obtain  the  necessary  data  and  the  inevitable 
wire  breakage*  forced  the  use  of  several  wires  during  the  measurements. 

All  data  were  obtained  from  traverses  along  the  Z-axis. 


There  are  two  types  of  analyses  of  the  hot-wire  output  which  are  of  interest. 
By  the  modal  analysis,  we  seek  to  break  down  the  fluctuations  received 
within  a  certain  bandwidth  into  contributions  from  each  flow  variable 
(modo);  e.g.,  the  velocity  or  the  density.  The  bandwidth  chosen  may  be 
small,  in  which  case  the  spectrum  of  each  separate  mode  can  be  obtained,  or 
large  enough  to  cover  the  entire  frequency  range,  in  which  case  the  fre¬ 
quency-integrated  contribution  of  each  mode  to  the  total  hot-wire  output 
can  be  determined.  In  the  former  instance,  the  spectral  analysis  of  the 
signal  yields  the  relative  importance  of  each  mode  at  different  frequencies 
(or  wave  numbers)  and  allows  for  the  computation  of  the  mode  autocorrelation 
function  and  its  turbulence  scales.  In  any  case,  the  first  basic  step  con¬ 
sists  of  measuring  the  total  wire  output  at  each  of  at  least  three  heating 
currents;  knowledge  of  the  mean-flow  properties  is  also  necessary  to  compute 
the  sensitivity  coefficients  (Reference  14).  Three  equations  are  thus 
obtained  and  one  can  solve  for  three  unknowns:  the  mass-flux  fluctuation, 
A(pu)/pu,  total  temperature  fluctuation,  AT/T,,  and  the  cross-correlation. 


"'Some  0.00005-inch-diameter  wires  survived  more  than  50  hours  of  tunnel 
time;  the  average  was  closer  to  20  hours. 
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r  of  these  two  fluctuations.  The  accuracy  of  this  measurement  can  be 
mt 

increased  if  more  than  3  measurements  are  taken  and,  in  fact,  15  different 
values  of  heating  current  have  been  chosen  in  the  present  experiment.  The 
highly  redundant  system  of  15  equations  was  then  solved  by  the  least-squares 
method . 

At  each  axial  station,  the  wire  was  first  placed  just  outside  the  wake,  its 
current  was  fixed  to  the  desired  value,  its  time  constant  was  measured  by 
the  square-wave  method,  and  the  compensating  amplifier  was  adjusted  to  that 
value  of  the  time  constant.  The  wire  was  then  traversed  radially  across 
the  wake  and  a  continuous  trace  of  its  average  voltage  and  the  frequency- 
integrated  mean  square  of  its  fluctuating  output,  deriving  from  the  maximum 
possible  bandwidth  of  the  system,  were  obtained.  This  procedure  was 
repeated  for  each  of  15  different  wire  currents  for  a  total  of  30  data 
curves  at  each  axial  station.  There  were  15  mean-wire  voltage  versus 
radius,  such  as  shown  in  Figure  6,  and  15  traces  of  the  frequency- integra ted 
mean-square  ac  component  of  the  wire  output  versus  radius.  An  example  of 
the  latter  is  shewn  in  Figure  7. 

If  the  frequency  response  of  the  wire-amplifier-compensator  system  were 
linear  over  the  range  of  the  flow  fluctuation  spectrum,  the  data  mentioned 
in  the  above  paragraph  would  be  sufficient  to  resolve  the  signal  into  the 
frequency-integrated  contributions  from  each  turbulence  mode  at  any  desired 
point  in  the  wake.  However,  neither  the  amplifier  nor  the  compensator 
could  attain  a  linear  transfer  function  over  the  500-kilocycle-per-second 
(keps)  range  where  signals  were,  by  preliminary  measurement,  known  to  occur. 
In  addition,  because  of  the  local  Reynolds  number  variation  in  the  wake, 
the  constant-current  wire  underwent  a  change  of  its  time  constant  with 
radius  and  thus  a  change  of  its  own  transfer  function  (frequency  response). 
Because  of  the  intense  turbulence,  proper  compensation  of  the  wire  within 
the  wake  was  not  possible.  Even  if  it  were,  the  point-by-point  adjustment 
of  the  compensator  transfer  function  would  increase  the  already  taxing 


16 


therefore  clear  that  s,orae  of  the  signals  contained  in  the  turbulent  fluid, 
(i.e.,  those  at  the  higher  frequencies),  were  missing  from  the  15  traverses 
previously  mentioned. 

To  overcome  this  difficulty,  the  following  principle  was  invoked:  the 
amount  by  which  an  imperfect  mechanical  system  distorts  (e.g.,  attenuates) 
an  input  signal  can  be  found  by  measuring  the  transfer  function  (i.e.,  the 
frequency  response)  of  the  system  itself  and  of  its  output.  Thus,  the  out¬ 
put  spectra  of  the  wire  at  each  axial  position  were  recorded  for  a  nunoer 
of  currents  and  radial  positions*.  The  results  were  combined  wi'h  tne 

known  overall  transfer  function  of  the  wire-amplifier-compensator  system 

2 

to  yield  the  ideal  frequency-integrated  mean  square  voltage  output,  e*  , 
of  the  wire  as  a  function  of  axial  and  radial  coordinates  x  and  y 


9  2  2  2  2 

where  e  is  the  recorded  mean  square  output,  W^G^/C  G  the  overa1!  transfer 

function,  V^(f)  the  wave  analyzer  output  at  each  frequency,  and  f  and  f^ 
are  the  1-kilocycle-per-second  analyzer  bandwidth.  Hie  error  ratio,  J, 
thus  obtained  was  then  applied  to  the  mean  square  output  from  the  constant- 
current  traverses  of  the  type  shown  in  Figure  7  to  correct  for  the  nonlinear 
response  of  the  system.  A  separate  computer  program  was  set  up  to  bring 
about  this  response  restoration  effect;  this  is  the  WEB-II  program  (See 
Appendix  B) .  Figure  9  shows  the  variation  of  the  wire  time  constant  and 
error  ratio  across  the  wake  17  diameters  downstream  of  the  base.  Interest¬ 
ingly,  J  approaches  a  limit  independent  of  wire  current  as  the  latter 


*This  was  done  for  each  of  8  currents  and  at  10  radial  positions  for  a 
total  of  80  spectra  at  each  axial  station.  Combined  with  the  mean  and 
turbulence  measurements,  a  total  of  115  traces  were  necessary  at  each 
axial  station.  An  example  of  such  spectra  is  shown  in  Figure  8. 
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FIGURE  9,  MEASURED  VARIATION  OF  WIRE  TIME  CONSTANT  (LEFT)  AND  ERROR 
RATIO  (RIGHT)  ACROSS  THE  WAKE  FOR  CURRENTS  VARYING  FROM 
6.425  MD  TO  2.976  MD.  CONDITIONS  ARE  X/D  =  17,  WIRE  3-4/2 
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increases.*  It  should  be  noted  that  the  time  constant,  M,  inside  the  wake 
was  computed,  for  each  current,  by 

«  '  Mi  [o4t 7 

where  is  the  time  constant  measured  at  the  wake  edge,  I  is  the  current, 
and  R  is  the  wire  resistance.  In  this  way,  reliance  on  a  theoretical  cal¬ 
culation  of  M  is  largely  circumvented  and  only  the  dependence  of  M  on 
overheat  and  Reynolds  number  is  retained  (see  Appendix  A).  A  final  point 
concerning  the  role  of  the  error  ratio,  i.e.,  the  need  for  the  response 
restoration  procedure,  is  shown  in  Figure  10.  We  see  that,  in  the  region 
investigated,  J  decreases  towards  the  far  wake  but  remains  higher  than 
unity.  Corrections  were  thus  necessary  even  at  about  100  diameters  down¬ 
stream. 

The  procedure  described  above  was  necessary  to  give  correct  results  for 
the  modal  analysis  of  the  data,  that  is,  correct  values  of  the  frequency- 
integrated  but  modally  resolved  fluctuations  from  point  to  point.  The 
computer  program  set  up  to  handle  the  above  "response  restoration"  of  the 
anemometer  system  (see  Appendix  B)  also  supplies  the  fluctuation  spectrum 
(i.e.,  the  true  distribution  of  the  flow  disturbances  with  frequency 
corrected  for  any  shortcomings  of  the  instrumentation)  at  each  axial  posi¬ 
tion  for  each  of  10  radial  locations  and,  at  each  of  the  latter,  at  each 
of  8  heating  currents.  The  output  of  spectra  listed  in  WEB-II 
(Appendix  B) ,  thus  formed  the  basis  of  the  spectral  analysis  wherein  the 
spectral  density  of  the  temperature  (density)  and  axial  velocity  fluctuation 
and  of  the  cross-correlation  could  be  resolved  at  each  point  in  the  wake. 


Figure  9  shows  that  the  zero-current  value  of  the  time  constant  in  the 
free-stream  approaches  the  theoretical  value  of  about  0.1  millisecond  by 
which  the  frequency  response  to  800  kilocycles  per  second  was  originally 
estimated.  Because  of  the  large  increase  of  the  time  constant  in  the 
wake,  however,  the  response  capability  decreases  there  and  results  in 
the  large  J's  measured. 
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FIGURE  10.  ERROR-RATIO  VARIATION  IN  THE  AXIAL  DIRECTION  AT  THE 
HIGHEST  WIRE  CURRENT  USED. 
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The  essential  point  hare,  hitherto  not  appreciated,  is  that  at  the  present 
state  of  the  art,  spectral  measurements  are  mandatory  whether  or  not  a 
spectral  analysis  of  the  data  is  intended. 

2.6  DATA  REDUCTION  AND  PROCESSING 
2.6.1  SIGNAL  PROCESSING 

It  is  now  necessary  to  describe  how  the  data  as  obtained  above  were  pro¬ 
cessed  through  the  electronic  equipment  shown  in  Figure  5.  This  descrip¬ 
tion  utilizes  the  principles  rather  than  the  exact  characteristics  of  the 
electronic  components  used  in  this  experiment  (the  latter  being  given  in 
Appendix  A)  so  as  to  be  of  value  to  others  engaged  In  similar  work. 

For  a  fixed  current,  the  sensible  real  input  to  the  measuring  system  of  a 
constant-current  anemometer  is  a  voltage  fluctuation  into  the  wire. 

Between  the  latter  and  the  recording  device,  the  signal  proceeds  through 
the  elements  of  the  block  diagram  of  Figure  5  and  is  thus  affected  by  the 
transfer  functions  of  these  elements.  Generally,  these  transfer  functions 
execute  selective  filtering  or  amplification,  squaring,  averaging,  inte¬ 
grating,  and  so  on.  Because  these  functions  are  often  unobtainable  from 
the  manufacturer,  they  were  measured  carefully  prior  to  this  work,  a  pro¬ 
cedure  necessary  to  increase  accuracy  and,  as  shown  in  the  preceding 
section,  to  make  the  wirk  possible  in  the  first  place. 

Consider  that  the  continuous  Fourier  spectrum  of  the  flow  into  the  wire 

(or,  otherwise  stated,  the  output  of  the  wire  if  the  latter  had  no  thermal 

lag)  is  the  sum  of  contributions  of  a  large  number  of  discrete  oscillators, 

each  vibrating  at  frequency  f,  and  contributing  a  mean  square  voltage 
2  2 

e  (f)  Flow  (volts)  /cps  in  the  interval  Af.  Then,  since  these  contributions 
are  additive  in  the  squete,  the  overall  wire  signal  is 


t 


Because  of  the  well-known  thermal  lag  phenomenon  (Reference  11),  the 

2  2 

actual  mean  square  wire  output  e  (f)  .  will  be  lower  than  e  (f),,  ,  and 

wire  v  ''flow’ 

the  transfer  function  W(f)  of  the  wire  is  then  defined  by 


W(f) 


e(f) 

e(f) 


flow 

wire 


£  1 


(4) 


The  signal  is  now  inserted  in  an  ac  amplifier  with  transfer  function  G(f) 

which  can  be  split  into  a  frequency-independent  zero-frequency  gain,  G  , 

and  the  amplifier  drop-off  characteristics  G/GQ(f);  the  amplifier  output, 

e(f)  can  then  be  defined  by 

amp ,  J 


G(f) 


e(f  ) 


amp 


e(f  ) 


wire 


G 

G 


(f)  G  <  G 
o  —  o 


(5) 


The  signal  next  proceeds  to  the  compensating  amplifier  which  has  usually  a 

zero-frequency  gain  of  unity,  a  transfer  function,  C(f),  and  an  output, 

e(f)  ,  such  that 
comp 


C(f) 


e(f) 


comp 


e(f) 


amp 


(6) 


Note  that  C(f)  is  also  a  strong  function  of  the  compensator  time-constant 
setting  and  that  it  generally  rises  well  above  unity  to  some  maximum  and 
thereafter  decays. 


If  Equations  (4),  (5),  and  (6)  are  combined,  the  mean  square  voltage  com 

2 

pensator  output  for  a  specific  flow  input,  e  in  the  passband  Af 

around  the  frequency,  f,  is  given  by 


e2(f)rnnn  Af 
comp 


C2CO  G2(f) 

w2(f) 


e2(f) 


flow 


Af  (volts)' 


(7) 


It  should  be  noted  that  the  quantity  CG  is  a  given  function  of  frequency 
for  a  given  setting  of  the  electronics.  The  quantity  CG/W,  on  the  other 
hand,  which  ;  l  equivalent  to  the  overall  transfer  function  (OTF)  is  addi¬ 
tionally  a  function  of  the  wire  time  constant  and  thus  a  function  of  the 
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wire  and  flow  characteristics,  changing  from  current  to  current  and  from 

point  to  point  in  the  wake;  however,  it  is  insensitive  to  the  magnitude  of 
2 

e  (f)fiow  within  broad  limits.  The  OTF  was  provided,  as  a  function  of 
frequency  at  each  wire  current  and  point  in  the  wake,  by  the  output  of  the 
WEB-11  program,  and  a  typical  plot  is  shown  in  Figure  11.  Note  that  the 
system  response  deteriorates,  although  not  greatly,  as  the  frequency 
increases. 


There  are  two  alternatives  to  further  processing  of  the  data.  For  spectral 
analysis,  one  uses  a  wave  analyzer  with  constant-bandwidth,  fg  cps;  in  this 
case,  the  analyzer  also  incorporated  a  detector  with  transfer  function,  H, 
such  that 


v2(f) 


and  thus 


H2  je2(f)  fR  (volts)4"  (dc) 

I  c  omp  d 


(8) 


e2(f) 


W 


flow 


2  2  2 
C  G  H  f 


B 


(volts)  (mean  square) 


(9) 


The  integrated  fluctuations  obtain  when  a  wideband  rms  voltmeter  obtains 
the  sum  of  the  compensator  output 


VTVM 

and  thus 

2 

6  VTVM 


i 


2  2 
e  (f)  df  (volts)  (mean  square) 
comp 


(10) 


Jo  W 2 


2  2 

e2(f).,  df  (volts)2  (mean  square)  (11) 
l  t  low 


This  voltmeter  also  has  a  detector  of  transfer  function  A  which  converts 
the  mean  square  fluctuating  voltage  into  a  dc  voltage  E'  such  that 


E'  "  A  eVTVM  V°ltS  (dc) 


(12) 


for  use  with  the  X-Y  recorder. 
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FIGURE  11.  TYPICAL  VARIATION  WITH  FREQUENCY  OF  THE  OVERALL  TRANSFER  FUNCTIONS 
AT  A  POINT  IN  THE  WAKE.  NOTE  IMPROVEMENT  IN  RESPONSE  AS  THE 
WIRE  CURRENT  I  IS  LOWERED. 


When  both  the  spectrum  and  the  integrated  measurements  are  performed,  we 
have  available  a  cross-check  of  the  experiment  since,  by  relating 
Equations  (9),  (11),  and  (12)  we  can  write 


EJ_  _  2 

A  *  eVTVM 


f  C2  G2 

Jo  W2 


e2(f) 


flow 


df  = 


This  check  was  performed  at  each  measuring  point 
program  WEB-II  and  gave  satisfactory  results,  as 

Returning  now  to  Equation  (11),  we  rewrite  it  as 


with  the  aid  of  computer 
can  be  seen  in  Figure  12. 


(14) 


For  an  ideal  amplifier  (flat  to  infinite  frequency),  G/Gq  is  unity;  C  is 
supposed  to  have  a  characteristic  exactly  inverse  to  W  below  some  frequency 
at  which  it  peaks,  and  if  e2(f)flow  -  0  beyond  that  point,  then  C/W  is 
also  unity  over  the  spectrum.  In  this  case,  the  ideal  mean  square  voltage 


VTVM 


I 


e2(f) 


flow 


df  (volts) ‘ 


(15) 


is  obtained,  and  thus  the  error  ratio,  J,  is  computed  as  explained  in 
paragraph  2.5. 


(16) 


Note  that,  as  explained  previously,  J  *  1  for  an  ideal  electromechanical 

system  (W/C  »  1,  G/Gq  =  1),  and  also  that  J  can  be  obtained  accurately 

if  H  and  f_  are  unknown  but  still  independent  of  frequency. 
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TYPICAL  AGREEMENT  BETWEEN  MODALLY  UNRESOLVED  RMS 
SQUARE  OUTPUT  AS  OBTAINED  BY  A  VOLTMETER  (VTVM2) 
AND  AS  RECONSTRUCTED  FROM  THE  SPECTRUM  (E02) ,  FOR 
DIFFERENT  CURRENT** . 
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A  last  point  concerns  the  subtraction  of  electronic  noise  from  the  measure* 
ment.  When  the  wire  heating  current  is  turned  off,  the  rms  voltmeter  reads 
a  noise  level  noise,  whereas  with  the  current  turned  back  on,  a 

reading  (e^^^)  meas  is  obtained.  Because  the  turbule'ce  and  electronic 
noise  are  uncorre la ted ,  one  obtains 


'VTVM 


(2  \  / 2  \ 

l  VTVM  I  meas  l 6 VTVM J  noi 


se 


(17) 


for  use  with  the  previous  formulas.  A  similar  correction  is  necessary  in 
the  spectral  measurement.  If  Vm  (f)  is  the  analyzer  detector  output  (in 
volts)  and  V^  (f)  is  the  corresponding  noise  at  frequency  f,  then 


is  the  corrected  reading  (in  volts  squared)  at  that  frequency. 


(18) 


Knowledge  of  the  wire  Reynolds  number  Re0  cause  from  the  diagnosis  of  the 
mean  properties  of  the  wave,  described  in  Reference  12.  The  modal  analysis 
into  the  rms  normalized  mass-flux  and  total  temperature  fluctuations  and 
their  correlation 


A(Pu)  (25) 

Pu 


rms 


(26) 


was  then  made  by  the  simultaneous  solution 
tions  of  the  type* 


AT 


rms 


(27) 


for  the  3  unknowns  of  15  equa- 


G2  e2 
o 


9  2  2  2. 

e  m  +  e  t  -  2  mtr  e  e 
m  t  mt  m  t 


(28) 


where  (e*2/G2  e2)  is  the  normalized  voltage  fluctuation  adjusted  for  the 
zero- frequency  amplifier  gain,  G0.  (The  bandwidth  of  these  modal  measure 
ments  was  in  each  case  sufficient  to  include  sensible  signals  at  all 
frequencies.)  Here  e*  (volts)  refers  to  the  ideal  frequency -Integra ted 
mean-square  wire  output  at  each  particular  value  of  the  heating  current 


“of^flov  df  <VOlt,)' 


(29) 


Note  the  similarity  with  Equations  (3)  and  (15),  which  say  that 

e*2  S  e**  =  E2  G2 

VTVM  o 


(30) 


*  Following  past  practice,  the  unknowns  are  obtained  by  the  computer  from 
the  shape  parameters  of  the  second-degree  curve  fitted  to  the  15  experi¬ 
mental  points.  (See  Appendix  D.) 
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2.6.2  MOu.'L  ANALYSIS 

Once  supplied  with  a  method  for  alleviating  the  response  shortcomings  of 
the  system,  the  turbulence  data  (i.e.,  those  described  by  the  15  traverses 
mentioned  earlier)  could  be  reduced.  The  sensitivity  coefficients  em  and 
e  to  mass  fiux  and  total  temperature  fluctuations,  respectively  were 
computed  by  the  following  slightly  modified  form  of  Morkovin's 
(Reference  14)  formulas 


=  e , 


Re  (5Nu 
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w  Nu  ^Re 
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where  cr  accounts  for  the  circuit  impedance,  T  and  K  are  related  to  the 

I  I 

resistance-temperature  slope,  and  Aw  is  Morkovin's  overheat  parameter 


w 


l  i  Sr 


(21) 


associated  with  the  wire  resistance  R  change  with  current,  I,  which  could 
be  directly  obtained  from  the  traverses  of  Figure  6.  For  the  type  of 
hot  wires  employed  (Reference  12),  the  zero-current  Nusselt  Number  NuQ 
based  on  the  stagnation  conductivity  is  given  by 


Nuq 

whence 

<)Nu 


0.04186  Re  +  0.1803  VRe7  +  0.1078 


0.04186  + 


0.0902 


Similarly 


0.99  + 


0.05661 


Re 
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(22) 
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Alternatively,  a  more  revealing  and  practically  important  modal  split 
decomposes  the  signal  into  vorticity  (velocity),  entropy,  and  sound 
fluctuations  (Reference  14  and  15).  It  is  significant  that  both  Kovasznay 
(Reference  15)  and  Kistler  (Reference  17)  have  confirmed  the  relative 
unimportance  of  the  sound  mode  inside  compressible  turbulent  shear  flows, 
compared  to  the  'orticity  and  entropy  modes.  The  data  of  the  present 
experiment  were  thus  further  reduced  on  the  assumption  that  there  were  no 
sound  fluctuations  inside  the  turbulent  wake.  At  each  ooint  in  the  wake, 
another  15  equations  are  obtained  of  the  type 


~r=2  "  4  j2  +  ea  °2  +  2  eT  ea  r0T  Ta 

G  e 
o 


where  the  vorticity  fluctuation 
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am  the  entropy  fluctuation 
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The  coeffocients  eT  and  ea  of  wire  sensitivity  to  vorticity  and  entropy 
fluctuations,  respectively,  have  been  derived  by  Morkcvin  in  terms  of  the 
corresponding  coefficients  em  and  e^.  given  by  Equations  (19)  and  (20) 

eT  -  em  +  <XeT  (35> 

e_  *■  Be  -  e 

O  r  T  m 


*  Following  past  practice,  the  unknowns  are  obtained  by  the  computer  from 
the  shape  parameters  of  the  second-degree  curve  fitted  to  the  15  experi¬ 
mental  points.  (See  Appendix  D.) 
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where  a  and  (3  are  functions  of  the  specific  heat  ratio  1  and  tne  Mach 
number  M 

-1 


a  s  l  + 


rtl  M2 
2 
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Again,  the  curve-fitting  process  was  used  to  solve  the  redundant  system 
of  15  equations,  such  as  (31),  for  T,  0,  and  the  correlation  ra„ . 

It  is  of  interest  to  look  at  the  variation  of  wire  sensitivity  -.c  the 
different  modes  with  heating  current,  as  shown  in  Figure  13.  The  out¬ 
standing  feature  is  the  existence  of  values  of  the  heating  current  for 
which  the  sensitivity  eT  to  vorticity  vanishes,  as  was  mentioned  in 
Paragraph  2.5.  This  important  fact,  first  recognized  by  Kovasznay 
(Reference  15),  implies  that  in  the  absence  of  sound,  the  density  fluctu¬ 
ations  alone  can  be  measured  directly  and  quite  simply  if  one  knows  the 
wire  characteristics  and  local  flow  properties  precisely,  in  the  present 
instance,  this  circumstance  was  utilized  to  provide  very  satisfactory 
cross-checks  of  the  entropy  fluctuation  results.  It  was  also  used  for 
studying  the  internment  nature  of  the  density  profile  of  the  wake. 

2.6.3  SPECTRAL  ANALYSIS  V 

As  mentioned ,  the\aim  of  Spectral  analysis  (or.  -ore  precisely,  -odal- 
apectral  analysls)\ls  to  show  t\  spectral  distributions  of  each  mode  of 
fluctuation.  The  hVr  to  this  an\lysls  is  that  the  equation  of  modal 
resolution.  Equ.tlon\(31) ,  is  independent  <-t  bandwidth.  Because  the 
fluctuations  are  constant  with  frequency  f  to  good  approximation  within 
the  lAilocycle-per-s*ond  bandwidth  used,  we  can  write  the  modal  equation 


SO?  -  4  T2(f)  +  4  02(f)  +  2  eT  efl  rOT(f)0(f)T(f) 

GZ 
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(39) 
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FIGURE  13.  VARIATION  OF  THE  WIRE  SENSITIVITY  COEFFICIENTS  AND  OVERHEATING 
PARAMETER  Aw  WITH  WIRE  CURRENT  AT  A  TYPICAL  POINT  IN  THE  WAKE. 


jl2  _  2 

at  each  frequence  f,  where  e  (f)/e  is  the  nondimensional  mean-square 
output  per  kilocycle.  The  sensitivity  coefficients  retain  their  previous 
values . 

By  way  of  clarification,  note  that  when  the  sum  of  signals  in  all  adjoining 
passbands  is  taken  and  the  result  integrated  over  frequency,  Equation  (39) 
gives 
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consistent  with  the  integrating  interval  df  (e.g.,  kilocycles  per  second. 

If  the  quantities  T(f)/ Tana  0(f)/ 0  are  defined  as  the  velocity  (vorticity) 

and  density  (entropy)  spectral  densities,  the  spectra  are  normalized  to 

unity;  that  is 
00 

T2(  f)df  -  1  (45) 

o 

a2(f)df  =  1  (46) 

The  spectra  can  next  be  normalized  by  the  characteristic  frequency  of  the 
turbulent  flow,  which  .is  u/A^  where  u  is  the  local  flow  velocity  and 
A  ^  is  the  integral  scale  appropriate  to  the  quantity  i  (AE  if  the  entropy 
scale,  A  T  the  vorticity  scale)  which  is  derivable  from  the  Eulerian  time 
integral  scale  for  the  single  wire;  the  latter  is  connected  with  the 
spectrum  (Reference  18)  so  that 
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SECTION  3 

PRELIMINARY  OBSERVATIONS 


3.1  THE  EXTERNAL  TURBULENCE  FIELD 

The  amount  of  free-stream  (ambient)  turbulence  in  the  wind  tunnel  was 
measured  prior  to  the  study  of  the  fluctuations  in  the  wake.  The  stream 
turbulence  is  often  a  critical  factor  in  that  it  can  add  to  the  noise  of 
the  electronic  equipment  in  obscuring  the  turbulence  level  in  the  wa>te 
itself.  Initial  measurements  without  the  model  in  place  showed  rather 
large  fluctuations,  on  the  order  of  1  percent  in  the  mass  flow  and 
0.3  percent  in  the  total  temperature.  Measurements  performed  later  with 
the  model  in  place  end  with  the  computerized  data  reduction  process 
described  in  the  previous  paragraphs  disclosed  the  above  computations  to 
be  larger  than  the  actual  turbulence  level  by  a  factor  of  about  3.  At  no 
time  was  the  stream  turbulence  level  a  serious  obstacle  to  the  fluctuations 
in  the  wake,  including  the  points  of  observation  farthest  downstream. 

3.2  OBSERVATIONS  ON  THE  TRANSITION  PROCESS 

Although  this  work  was  not  aimed  at  investigating  the  events  in  the  laminar 
wake  leading  up  to  transition  to  turbulence,  some  interesting  phenomena 
were  detected  in  the  laminar  and  pretransitional  flow.  Briefly,  a  sharp 
energy  concentration  at  about  20  kilocycles  per  second  was  observed  to 
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occur  in  the  laminar  wake,  with  the  rest  of  the  spectrum  completely  devoid 
of  fluctuations.  As  the  wire  approached  the  transition  rone,  the  energy 
prominence  first  increased  in  magnitude  and  then  rapidly  began  to  feed 
energy  into  the  higher  frequencies.  Further  downstream,  the  peak  decayed 
so  that  a  completely  random  spectrum  remained  in  the  turbulent  wake.  This 
result  was  obtained  both  by  first  keeping  the  wire  stationary  and  moving 
the  transition  forward  toward  thn  probe  by  increasing  the  tunnel  pressure, 
and  also  by  keeping  the  flow  fixed  and  moving  the  wire  downstream  in  the 
laminar  wake  toward  the  transition  zone.  It  was  also  obtained  with  differ¬ 
ent  wires.  These  facts,  added  to  the  care  exercised  to  recognize  and 
reject  extraneous  signals,  strongly  suggest  that  the  phenomenon  is  real 
and  is  further  associated  with  the  transition  process.  The  intriguing 
possibility  of  finding  a  connection  between  laminar  instability  and  transi¬ 
tion  to  turbulence  suggests  that  further  exploration  of  this  phenomenon  is 
needed.  These  observations  are  pictured  in  Figure  14. 
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FIGURE  14.  UNPROCESSED  HOT-WIRE  SPECTRA  TAKEN  AT  A  FIXED  LOCATION  IN 
THE  WAKE  FOR  DIFFERENT  TUNNEL  STAGNATION  PRESSURES  SHOWING 
THE  ONSET,  INCREASE  AND  "SPECTRAL  SPREADING"  OF  L  AMINAR 
INSTABILITY. 
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SECTION  4 


SIMILARITY  CONSIDERATIONS 


We  shall  preface  the  description  of  the  experimental,  results  by  outlining 
current  empirical  predictions  of  the  distribution  of  turbulence  properties 
in  the  wake.  Our  chief  premise  is  that,  following  a  certain  relaxation 
period,  the  turbulence  properties  may  soon  adjust  to  a  self-preserving 
(similar)  form  connected  simply  with  the  local  mean  properties.  If  it 
exists,  this  similarity  behavior  should  he  closely  connected  with  the 
behavior  of  incompressible  turbulent  flows.  By  providing  experimental 
proof  of  a  more  general  rule,  insensitive  to  compressibility  effects,  we 
have  a  good  chance  of  extending  our  present  results  to  cover  the  dynamics 
of  wake  flows  at  very  high  Mach  numbers . 

The  dynamic  equilibrium  hypothesis  (DEH) ,  discoursed  at  length  by  Townsend 
(Reference  1) ,  prescribes  that  the  large  eddies  maintain  a  dynamical 
balance  between  the  extremities  of  turbulent  shear  zones  in  a  way  that 
simply  relates  the  fluctuations  and  the  mean  properties.  These  ideas, 
formulated  in  context  with  incompressible  flows,  have  been  recently 
extended  to  compressible  turbulent  flows  by  Morkovin  (Reference  10)  and 
Laufer  (Reference  11).  Morkovin  observed  that,  with  appropriate  scaling, 
the  turbulent  energy  distribution  in  sheai  flows  is  insensitive  to  geometry 


-39- 


and  Mach  number.  Extended  to  the  compressible  turbulent  wake  by  Laufer, 
and  restricted  to  adiabatic  flows,  this  observation  has  the  following 

implications : 

(1)  There  are  no  fluctuations  in  the  flow 
total  temperature 

(49) 

AT  =  o 
o 

(2)  In  consequence,  the  entropy  (density)  and 
vorticity  (velocity)  fluctuations  are  related 
through  the  local  Mach  number,  M 


(Ap/p  ) 


rms 


(Au/u) 


=  (7  -  1)  M 


rms 


(3)  The  temperature  and  velocity  fluctuations 
are  perfectly  anti-correlated 


rOT 


(4)  The  nondimens ional  turbulent  energy  is  a 

universal  function  g  of  a  properly  transformed 
transverse  coordinate  !)• 


(50) 


(51) 


(52) 


p(Au)Z  -  P'  u  scn> 

where  p'  and  U  are  appropriate  density  and 
velocity  scales,  respectively. 

r,:n  “r— ..  -  -»r 

.  Mnn  /52)  As  Laufer  (Reference  11)  pointed  out,  a  more- 
arising  from  Equation  (52) .  As  uauie  v 

general  form  of  Equation  (52)  is 


/a  \ 2  vm 

p(Au)  Y 


-  p'  U2  Ym  g(T)> 


(53) 


where  m  =  0  or  1  for  planar  or  axi- symmetric  geometries,  respectively,  and 
where  the  quantity  AQ  stands  for  the  root-mean- square  (dimensional)  fluc¬ 
tuation  of  the  quantity  Q.  Here  Y  is  a  transformed  radial  coordinate 
which  might  be  of  the  type 

-2  rV  P 

Y  =  2  f  f  Y  dr  (54) 

■'O  oO 

We  immediately  see  that  the  latter  choice  is  unacceptable  if  p  (Aj)2  is  to 
be  truly  similar,  because  we  would  then  have 

~^2~  =  |  (X  Tl)  g(T))  (55) 

p'  IT  Y 

because  p/p^  is  generally  a  function  of  both  X  and  Y.  We  therefore  choose 
Y  =  Y;  that  is,  we  use  the  premise  that  similarity  in  the  turbulent  energy 
exists  in  the  physical  rather  than  the  Howarth-Dorodnitzyn  plane.  This 
also  implies  a  certain  insensitivity  to  the  geometry,  in  agreement  with 
Morkovin's  observations.  Thus  Equation  (53)  reduces  to  Equation  (52). 

We  will  next  recognize  that  the  scales  are  given  by 

Density  Scale:  p'  =  p  (56) 

~  '00 

Velocity  Scale:  U  =  u  -  u(o)  (57) 

were  oo  refers  to  conditions  outside  the  wake  and  (o)  refers  to  conditions 
on  the  axis.  The  first  obvious  result  is  that 


_P  (Au)2 

u2 


gOl) 


which  for  the  axis  reduces  to 

.2 


P  (o)  (Au(o)) 


g(o) 


00 


u‘ 


■  constant 


(58) 


(59) 
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in  contrast  with  the  incompressible  wake  for  which  Au(o)/U  =  constant. 

Since  p  (o)  <  pw  in  the  compressible  wake,  the  fluctuations  are  initially 
higher  (and  much  more  so  for  hypersonic  wakes)  than  for  low-speed  wakes. 
Eventually,  of  course,  p(o)  ■  p '  and  the  constant  ratio  of  g(o)  is  attained. 

The  radial  variation  of  the  velocity  fluctuations  in  terms  of  the  variable 
1  |  (60) 


where  L  is  the  transverse  wake  scale  (to  be  discussed  shortly),  obtains, 
directly  from  Equations  (58)  and  (59) 


=  =  f(n) 

g(o) 


(Cl) 


where  f (T))  is  another  unique  function  of  7).  The  universality  of  this 
function  can  be  tested  by  comparing  (for  example)  velocity  fluctuation 
profiles  (Au)2/ (Au(o) )2  in  incompressible  flows  with  the  profiles 
p  (Au)2/p(o)  (Au(o))  of  the  present  experiment.  At  any  rate,  since 
p  >  p(o),  the  profile  (Au)/(Au  (o))  in  the  compressible  wake  is  somewhat 
thinner  than  in  the  incompressible  wake. 

The  following  corollaries  can  be  drawn  from  Equations  (58)  and  (61) 


_L _ £_  ^Au 

w2 


g  CO) 


Pu 


(Au/u)' 


p (o)  u(o) 2  (Au(o)/u(o))2 


f  (7?) 


where  w  i3  the  so-called  velocity  defect 


(62) 


(63) 


w 


u  -  u(o) 
00 


u, 


00 


and  also 


(64) 


Au(o)  2 
u(o) 


_ _ 

(1  -  r)  (1  -  w)2 


(65) 
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where  the  density  defect 


p  (o) 


As  will  be  discussed  below,  the  defects  W  and  r  are  known  functions  of  the 
axial  coordinate  X  so  that  the  X-dependence  of  the  velocity  fluctuations 
can  be  predicted.  In  the  meantime,  we  again  note  the  influence  of  com¬ 
pressibility  on  the  axial  decay  of  the  velocity  fluctuation  on  the  axis, 
as  given  by  Equation  (66)  and  the  axis  value  of  Equation  (62) 

/  \  ?  n  u  2 


Au(o) 

u„ 


2  p  v* 

e<°>”  T?)  "  «<°> 


In  both  cases,  we  observe  that  the  fluctuation  decay  is  faster  ir.  the  com¬ 
pressible  wake  because  of  the  decrease  in  r  with  axial  distance. 

4.2  TEMPERATURE  AND  DENSITY  FLUCTUATIONS 

The  prediction  of  Equation  (49)  for  adiabatic  flows  together  with  the 
likelihood  of  insignificant  sound  fields  allows  us  to  construct  the  simi¬ 
larity  behavior  of  the  temperature  and  density  fluctuations 

_  .  dT  (68) 

p  T 

iA£l  „  (69) 

p  T 

where  the  parentheses  again  denote  rms  values.  Equation  (50)  rewritten  as 

-  CY  -  1)  M2  —  (70> 

T  '  u 

immediately  leads  to  the  radial  profile  of  the  temperature  fluctuations  as 
normalized  to  the  axis  value 


(nisi) 

l  AT  ) 

P(o) 

\  u  ) 

\AT(o )/ 

and  its  corollaries 


ELsll 

(aL°l) 

P3 

\  U  / 

EL  q) 

P 

\  u  f 

P(o) 

(aLol) 

p 

\  U  / 

-m.  y  . 

(Ap>  (o)/ 

_ c  m/p 
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f(TJ) 


-  f  (▼)) 

f(T)) 


(72) 

(73) 

(74) 


Here  we  note  the  analogies  among  Equations  (61),  (71),  and  (72).  In  the 
same  spirit,  we  inquire  about  the  axial  variation  of  the  axis  values  of  the 
density  and  temperature  fluctuations.  In  analogy  to  Equation  (59)  we 
obtain 


&T(o) 

T(o)  -  Too 


g(o) 


(7  -  1)  M(o) 


2  «  (1  +  e 

e  (i  -  w) 


(75) 


where  Q  is  the  temperature  defect 


6  *  T-^°j  '  Tfia  (76) 

ioo 

and  where  M(o)  is  the  Mach  number  of  the  axis.  A  more  convenient  expres¬ 
sion  can  be  given  using  the  experimentally  determined  (Reference  12) 
relation 

0  a  (7-1)  M^w  S  (77) 

where  S  is  a  known  number  on  the  order  of  0.8  related  to  the  turbulent 
Prandtl  number.  Then 


ano) 

T(o)  -  Too 


1 

S 


g(o) 


pi  o) 


1/2 


(1  -  w) 


(78) 


Thu  ratio  of  the  temperature  fluctuations  on  the  axis  to  the  temperature 
difference  across  the  wake  is  therefore  not  constant,  as  Gibson  (Reference  4) 
found  for  the  sphere  wake  in  water,  but  losely  follows  the  behavior  of 


the  density  and  velocity  defects.  In  fact,  the  density  fluctuations 
given  by 

^  p(o)  -  i  [*<°>  -%f]  1/2<‘  -  U>  <79> 

analogous  to  Equation  (78),  obvious  y  increase  with  axial  distance. 

Both  Equations  (78)  and  (79)  also  predict  that  at  very  far  distances 

from  the  body  (p(o)  ■  p  ,  w  ■  o) ,  the  temperature  and  density  fluctuations 

00  1/2 

so  defined  approach  the  value  (g(o))  /S. 


The  behavior  of  the  fluctuations  as  given  above  carry  (concealed)  the 
surmise  that  the  appropriate  density  sca’.e  p'  ,  is,  in  fact,  given  by  p'  «  p 
(see  Equation  (56)).  Furthermore,  the  constant  g(o)  is  needed  in  applying 
these  formulas.  An  interesting  relation  devoid  of  these  requirements  is 
the  ratio 


/  AP(o)  \  /  (  AU(o)  \  .  (l-.rj(l-w) 

V  Poo  -  p(°)/  /  \Uoo  -  u(°)/  S 


(80) 


obtained  from  Equations  (59)  and  (79).  Both  sides  of  Equation  (80)  have 
been  measured  Independently  in  this  experiment  and  a  check  would  be  much 


in  order. 


4.3  SPATIAL  DISTRIBUTION 

The  velocity,  density,  and  temperature  fluctuations  are  given  above  in 
terms  of  the  mean-wake  properties  and  particularly  in  terms  of  the 
defects  w,  6,  and  r.  This  representation  is  powerful  in  that  it  bypasses, 
an  exact  knowledge  of  the  mean-property  distribution  in  the  wake. 
Nevertheless,  the  latter  distributions  are  now  well  understood  for 
compressible  wakes  in  general  and,  by  the  results  of  Reference  12,  for 
the  wake  at  hand  in  particular.  To  describe  the  functional  forms  AuCn,X) » 
etc.,  we  will  first  list  the  main  features  of  the  mean  flow  field  of  this 
wake  as  given  in  Reference  12. 


The  coordinates  a*"e  defined  by 


Axial  coordinate  X 


X  -  X 
_ o 

/V 


(81) 


Radial  coordinate  T)  =  “■  <82^ 

Here  X  is  the  actual  distance  behind  the  body  base,  Xq  the  X-position  of 
the  effective  wake  origin  (found  to  lie  about  8  diameters  behind  the  base), 
and  (CdA)^2  is  the  virtual  (drag)  diameter  of  the  wake;  in  the  present 


instance* 


<cda) 


1/2 


0.216  cm 


The  wake  transverse  scale,  L,  is  given  by 


L  - 


V 


1  1/2 


UttwJ 

and  the  velocity  defect 

2/3 


w  = 


Tfl/3 


(¥) 


■=2/3  ’ 


12.8 


(83) 


(84) 


It  should  be  noted  that  the  T]  of  Reference  12  was  formed  by  using  the 
Howarth-Dorodnitzyn  radius,  which  was  very  similar  to  the  physical  radius 
Y,  beginning  a  very  short  distance  away  from  the  axis. 

The  density  and  temperature  defects  are 

0  -  (7  *  1)  ««  «  S  (8 


*The  value  (CdA)1/2  »  0.278  cm  given  in  Reference  12  is  an  average 
Including  measurements  in  the  up-down  direction,  less  accurate  than  t  e 
value  0.216  cm.  However,  the  latter  value  was  used  to  compute  X  in 


Reference  12. 


and 
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e 

9  +  1 


(86) 


Again,  the  quantity  aT2/3  appearing  in  place  of  S  in  Equation  (38)  in 
Reference  12  has  been  replaced  by  a  more  general  function,  S,  of  the  tur¬ 
bulent  Prandtl  number,  aT,  which  varies  from  about  0.78  to  0.88  along  the 
wake,  tending  toward  0.84  far  along  the  wake  (see  Figure  35,  Reference  i2). 


The  radial 

u  ~ 


T 


P  15 


variation  of  properties  is  as  follows 


(87) 


(88) 


(89) 


The  values  of  w,  Q,  and  r  found  experimentally  in  Reference  12  are  listed 
in  Table  I.  It  will  be  noted  that  they  obey  the  above  relations 
(Equation  (64),  for  example)  only  after  a  certain  relaxation  length  behind 
the  transition  region.  We  therefore  anticipate  that  a  similar  relaxation 
oerlod  will  be  found  in  the  fluctuations,  i.e.,  that  the  findings  of 
Paragraphs  4.1.1  and  4.1.2  may,  if  at  all,  fit  the  data  only  after  some 
distance  behind  the  body.  Based  on  Reference  12,  this  distance  was  a 
minimum  of  40  in  X,  contrasted  with  the  distance  X  at  90  covered  by  the 
measurements. 

By  combining  the  equations  of  Paragraphs  4.1.1  and  4.1.2  with  those  of 
Paragraphs  4.1.3,  the  behavior  of  the  fluctuations  tlong  the  wake  can  be 
predicted.  For  brevity,  these  combinations  are  outlined  in  Table  II.  At 
Intermediate  X's  the  similarity  behavior  of  the  compressible  wake  differs 
distinctly  from  that  of  the  incompressible  wake  because  of  density  factors 
such  as  the  ratio  p  (o )!pV)  and  the  defect,  r.  This  intermediate  region 


-47- 


TABLE  I 


MEAN-WAKE 

PROPERTIES 

on  wake  axis 

(FROM  REFERENCE 

12) 

X-Station 

X 

w 

e 

r 

0 

14.71 

0.276 

0.754 

0.433 

1 

20.60 

0.222 

0.622 

0.379 

2 

26.50 

0.181 

0.522 

0.334 

3 

32.39 

0.150 

0.435 

0.303 

4 

38.28 

0.137 

0.327 

0.253 

5 

44.17 

0.107 

0.313 

0.246 

6 

50.07 

0.0951 

0.283 

0.225 

7 

55.96 

0.0731 

0.231 

0.186 

8 

61.85 

0.0688 

0.211 

0.167 

9 

67.75 

0.0615 

0.181 

0.165 

10 

73.64 

0.0595 

0.171 

0.168 

11 

79.53 

0.0564 

0.173 

0.133 

12 

85.42 

0.0559 

0.168 

0.116 

13 

91.32 

0.0490 

0.149 

0.138 

should  not  be  confused  with  the  relaxation  region  mentioned  previously; 
we  do  still  expect  similarity  behavior  with  radial  distance  in  the  former 
but  not  in  the  latter.  Again,  note  that  Au/(uw  -  u(o))  and  &p/ (pm  -  p(o)) 
change  with  X  and  the  former  decreases  toward  its  asymptotic  limit  ^g(o) , 
whereas  the  latter  increases  toward  its  limit  ^g(o)/S--limits  which  will, 
incidentally,  be  identical  if  the  Prandtl  number  (and  thus  also  S)  is  unity. 
The  velocity  fluctuation  decay  obeys  incompressible  rules*  only  in  the 
limit  of  very  large  X. 

It  remains  to  evaluate  g(o)  which,  in  view  of  our  committment  to  the  dynamic 
equilibrium  hypothesis,  should  be  similar  to  the  incompressible  values. 
Unfortunately,  the  axi- symmetric  experiments  of  References  2,  3,  and  5 
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lack  the  adequacy  and  clarity  of  data  from  which  g(o),  which  seems  to  tend 
towards  unity  near  the  body,  can  be  obtained  at  large  X.  It  is  interesting, 
however,  that  Townsend  (Reference  1,  Page  135)  obtained  the  value  g2  (D)  = 

0.13S  in  the  two-dimensional  wake. 

We  close  this  discussion  with  the  remark  that  the  X  dependencies  of  the 
velocity  defect,  the  transverse  scale  L,  etc.,  as  shown  above,  have  the 
exact  form  exhibited  at  low  speeds,  including  the  numerical  values  of  the 
coefficients  (RT  for  example).  Furthermore,  these  expressions  have  been 
validated  in  the  compressible  wake  by  the  experimental  results  of  Reference  12 
and  thus  add  to  the  expectation  that  the  predictions  of  the  DEH  will  also  be 
validated  for  the  fluctuation. 

4.4  CORRELATIONS,  SPECTRA,  AND  SCALES 

There  are  some  definite,  although  semiqja5 itative ,  predictions  the  DEH  makes 
concerning  the  local  velocity- temperecure  correlation  (cross-correlation) 
as  well  as  the  fluctuation  spectral  density.  If  our  picture  of  large-eddy 
equilibration  is  correct,  then  we  expect  Equation  (51)  to  hold,  and  princi¬ 
pally  so  for  the  large  eddies  (low  frequencies);  the  rationale  is  that  large 
eddies,  unaffected  as  they  are  by  internal  mixing,  best  represent  the 
extreme  conditions  across  shear  zones.  In  fact,  velocity- temperature  cor¬ 
relation  measurements  (References  10,  15,  and  17)  show  that  these  two 
quantities  are  usually  strongly  (although  not  perfectly)  anti-correlated. 

By  the  same  token,  one  expects,  in  the  absence  of  £T0  fluctuations, 

Equation  (50)  for  the  temperature-velocity  fluctuation  ratio  to  perform 
best  at  the  lower  frequencies.  Both  these  phenomena  are  illustrated  by 
Morkovin  (Reference  10)  for  the  boundary  layer. 

The  Integral  scale  of  turbulence  (i.e.,  the  autocorrelation  macroscale) 
should  be  expected  to  scale  with  the  transverse  scale  of  the  wake  mean  flow 

A  **  L 

with  the  constant  of  proportionality  to  be  determined  by  the  experiment. 

The  microscales  of  the  autocorrelation  functions  will  be  discussed  in  a 
later  report. 
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SECTION  5 


RESULTS  OF  MODAL  ANALYSIS 


The  experimental  results  are  summarized  In  the  form  of  the  axis  values  in 
Table  III.  In  this  section  we  shall  present  the  results  of  the  modal 
analysis  of  the  wake  turbulence  (the  frequency- Integrated  fluctuation 
Intensities)  beginning  with  the  wire-sensible  variables,  l.e.,  the  fluctu¬ 
ations  In  the  mass-flow-flux  pu  and  the  total  temperature  T0,  along  with 
their  cross-correlation  coefficient.  These  results,  put  out  by  the 
WEB-IV  program  (Appendix  D)  derive  from  the  solution  of  Equation  (28) 
(Paragraph  2.6.2)  with  the  inputs  of  the  latter  deriving  from  the  experiment. 

Figure  15  shows  the  variation  of  the  rms  mass-flux  and  total  temperature 
fluctuations  along  the  axis,  where  the  fluctuations  are  normalized  by  the 
local  (axis)  mean  values.  Note  that  the  total  temperature  fluctuations  are 
in  fact  quite  small,  attaining  a  maximum  value  of  7  percent  In  the  region 
studied  and  decaying  to  about  1  percent  at  100  diameters  downstream  of 
transition.  This  is  considerably  below  all  other  fluctuation  levels  and 
especially,  as  will  be  seen  below,  the  fluctuations  In  the  static  temper¬ 
ature.  We  already  see  an  Indication  that  Morkovin's  strong  Reynolds 
analogy  (Reference  10)  for  adiabatic  flows,  which  implies  the  criterion  of 

Equation  (49),  is  at  work.  The  correlation  coefficient  r  between  mass  flux 
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FIGUU  16.  TYPICAL  RADIAL  VARIATION  OF  TUI 
WITH  THI  LOCAL  MASS  FLUX. 
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TABLE  III 

TURBULENCE  PROPERTIES  ON  WAKE  AXIS 


X-Statlon  X 

Au/u 

A  p/  p  A 

u/(ul<ru(0)) 

Ap/(poo-p(°)) 

0 

14.71 

0.210 

0.142 

0.549 

0.186 

1 

20.60 

0.156 

0.148 

0.548 

0.243 

2 

26.50 

0.113 

0.136 

0.511 

0.27? 

3 

32.39 

0.0937 

0.154 

0.532 

0.356 

4 

38.28 

0.0793 

0.125 

0.501 

0.369 

5 

44.17 

0.0531 

0.114 

0.442 

0.350 

6 

50.07 

0.0230 

0.100 

0.220 

0.344 

7 

55.96 

0.0361 

0.0993 

0.458 

0.434 

8 

61.85 

0.0350 

0.0852 

0.473 

0.424 

9 

67.75 

0.0260 

0.080 

0.430 

0.420 

10 

73.64 

0.0252 

0.0781 

0.399 

0.389 

11 

79.53 

- 

0.0649 

- 

0.423 

12 

85.42 

0.0224 

0.0658 

0.378 

0.500 

13 

91.32 

0.0192 

0.0591 

0.373 

0.368 

and  total 

temperature. 

also  shown  in  Figure 

15,  implies 

strong  correlation 

near  the  transition  zone,  but  the  fluctuations  become  totally  uncorrelated 
further  downstream.* 


The  variation  of  Apu/pu  and  ATq/To  with  radius  Y  is  shown  in  Figures  16  and 
17,  typically  represented  by  a  radial  survey  about  the  middle  of  the  wake. 


The  modal  analysis  according  to  Equation  (31)  has  yielded  the  local  tempera¬ 
ture  (density)  and  axial  velocity  fluctuations  normalised  to  their  respective 
mean  values,  the  same  quantities  normalized  to  the  local  mean  differences 


*Slnce  cross-correlations  should  be  ±1  in  unsteady  laminar  flows,  it  is 
suggested  that  the  variation  of  r  in  0  X  40  is  indicative  of  the  pro¬ 
gression  from  the  transitional  to  the  fully  established  turbulent  regime; 
this  agrees  with  the  relaxation  (X» 40)  distance  of  Reference  12.  However, 
this  is  not  to  be  confused  with  the  expected  behavior  of  r0T>  about  which 
will  be  said  later. 
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FIGURE  17.  TYPICAL  RADIAL  VAR  I AI I  OH  OF  THE  RMS  TOTAL  TEMPERATURE 
FLUCTUATION  NORMALIZED  WITH  THE  LOCAL  TOTAL  TEMPERATURE, 
COMPARE  MAGHITUDRS  WITH  FIGURE  16. 
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across  the  wake,  and  the  ve loclty-temperature  correlation  coefficient  which 
is,  of  course,  equal  in  magnitude  but  opposite  in  sign  to  the  velocity- 
dens  '  ty  correlation.  These  quantities  are  available  at  the  computer  output 
plotted  versus  the  physical  radius  Y  in  inches  as  well  as  the  normalized 
radial  distance  T)  -  Y/L  where  L  is  the  transverse  wake  scale  available  from 
the  mean  measurements  of  the  WEB-V  computer  output  (Reference  12). 

These  turbulence  data  will  be  discussed  first,  as  plotted,  without  regard  to 
the  DEH  and  its  extension  to  compressible  flows.  We  note,  first,  the  repre¬ 
sentative  plots  of  Figures  18,  19,  20,  and  21.  Both  modes  peak  at  some 
distance  off  the  axis,  with  the  peak  in  the  density  fluctuations  higher  and 
farther  out  from  the  axis.  In  terms  of  their  local  values,  the  velocity 
fluctuations  are  considerably  lower  than  the  density  fluctuations;  since  the 
Mach  number  in  the  wake  ranges  between  1.5  and  3  and  the  total  temperature 
fluctuations  are  very  small,  this  behavior  Immediately  confirms  the  general 
trend  of  Equation  (49).  In  terms  of  the  differences  across  the  wake,  these 
two  fluctuations  are  very  similar  and  their  maximum  values  are  about  50  per¬ 
cent  of  these  differences,  in  contrast  with  the  value  of  100  percent  which 
one  would  be  led  to  believe  from  rough  dimensional  arguments.  The  decay  of 
the  axis  velocity  and  density  fluctuations  (normalized  with  their  respective 
local  mean  values)  with  axial  distance  is  shown  in  Figure  22. 

We  can  best  illustrate  the  compressible-incompressible  dissimilarities  by 
summarizing  the  above  discussion  by  the  plots  shown  in  Figures  23  and  24. 

In  the  former,  the  rms  axial  velocity  fluctuation  normalized  with  the  exter¬ 
nal  velocity  is  shown  to  decrease  along  the  axis.  As  predicted  in  Table  II, 
this  fluctuation  decays  faster  than  the  (-2/3)  power  of  X;  this  should  be 
compared  with  the  incompressible  behavior  (References  2,  3,  and  5).  At  this 
stage  the  cause  can  be  attributed  to  the  density  factor  (1  -  r)  (cf.  Table  II 
and  Equation  (67))  which  of  course  Increases  with  X.  In  Figure  24,  the  rms 
density  and  axial  velocity  fluctuation  values  on  the  axis  are  shown  as 
functions  of  X,  normalized  to  their  respective  differences  pw  -  p (o)  and 
um  -  u(o)  across  the  wake.  Despite  the  data  scatter,  it  is  quite 
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FIGURE  24.  AXIAL  VARIATION  OF  THE  IMS  AXIS  VELOCITY  (TOP)  AND  DENSITY  (BOTTOM) 

FLUCTUATIONS  NORMALIZED  WITH  THEIR  RESPECTIVE  DIFFERENCES  ACROSS  THE  WAKE 


t  • 


obvious  that  these  quantities  are  each  far  from  unity  and  in  fact  change  in 
opposite  directions--the  velocity  fluctuations  decreasing,  and  the  density 
fluctuations  increasing,  with  increasing  X.  A  quick  look  at  Table  II 
indicates  that  this  behavior  agrees  well  with  the  dynamic  equilibrium 
hypothesis . 


Still,  the  central  test  is  the  behavior  of  the  rms  fluctuations  normalized 
to  the  axis  fluctuation  with  radial  distance;  these  are  plotted  against  the 
normalized  radius  in  Figures  25  and  26  for  X  >44.  For  the  velocity 
fluctuations,  the  apparent  lack  of  perfect  similarity  is  due  largely  tc 
experimental  scatter.  The  role  of  similarity  behavior  is  shown  better  in 

9 

the  density  fluctuations  (Figure  26)  and  here  it  will  be  seen  that  by 
replotting  these  points  according  to  Equation  (72),  their  dispersion  is 
considerably  decreased.  This  dispersion  is  much  worse  for  X  <  44  than 
Figure  26  indicates;  the  peal-  n  each  curve  of  Ap/Ap(o)  exceeds  2.0  nearer 
the  body.  A  further  point  worth  retaining  from  Figures  25  and  26  is  that 
the  peak  in  the  density  .fluctuation  is  farther  out  (near  T)  ■  2)  than  it  is 
for  the  velocity  (cf.  Figure  25).  This  difference  is  too  large  to  be 
explained  by  arguments  of  a  "thicker  temperature  wake"  and,  as  will  be  seen 
shortly,  practically  vanishes  when  the  data  are  replotted  by  the  method 
suggested  by  the  DEH. 
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A  test  of  Equation  (50)  which  enunciates  Morkovln's  strong  Reynolds 
analogy  (Reference  10)  (a  corrollary  of  the  DEH)  was  also  made  by  plotting 
the  velocity  temperature  correlation,  r^.  In  Figure  27.  The  r^versus  t) 
plot  presents  about  250  experimental  points  and  shows  that,  with  very  few 
exceptions,  r^  ranges  from  >0.8  to  -1.0.  No  systematic  radial  variation 
was  found,  but  anti-correlation  weakened  somewhat  (from  -1.0  to  -0.8)  with 
Increasing  X.  This  behavior,  already  hinted  by  the  weakness  of  the  total 
temperature  fluctuations  ATq,  Indicates  that  the  strong  Reynolds  analogy 
was  very  active  in  the  wake.  Note  tha^  Kovasznay  (Reference  15)  and 
Klstler  (Reference  17)  both  found  an  r^  of  about  -0.7  almost  constant 
with  distance  across  a  supersonic  turbulent  boundary  layer.  Physically 
Interpreted  anti-correlation  Implies  that  the  wire  Is  struck  alternately 
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FIGURE  25.  RADIAL  DISTRIBUTION  OF  THE  RMS  AXIAL  VELOCITY  FLUCTUATIONS  NORMALIZED  WITH  THE 
AXIS  VALUE,  FAR  IN  THE  WAKE. 
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FIGURE  26.  RADIAL  DISTRIBUTION  OF  THE  RMS  DENSITY  FLUCTUATIONS  NORMALIZED  WITH  THE  AXIS 
VALUE,  FAR  IN  THE  WAKE. 


FIGURE  27.  AXIAL  AND  RADIAL  VARIATIONS  OF  THE  VELOCITY -TEMPERATURE  CROSS-CORRELATION 
COEFFICIENT.  THE  MEAN  PROPERTY  VARIATIONS  ARE  SHOWN  AT  THE  TOP  FOR 
COMPARISON. 
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*  Old  and  fast  or  hot  and  slow  moving  fluld,  Illustrating  the  father 

ap  ranafet  across  the  sheer  rone  and  the  long  roach  of  the  large 
eddies.  ° 

SO  far,  the  data  have  been  purposefully  presented  and  discussed  in  a 
manner  showing  the  Inadequacy  of  low- speed  arguments  to  cover  the 
behavior  of  the  compressible  turbulent  wake.  If  these  arguments  are 
indeed  special  cases  of  the  DEH  in  the  constant-density  limit,  then 

d  ‘“r  e"8‘"eer1"8  COr"U“°"»  «  acbleved  by  regrouping  the  present 
data  aa  indicated  by  the  formulas  of  Paragraph  4.1.1.  shall  present 
these  correlations  below. 

First,  we  look  for  similarity  In  the  radial  profile,  of  velocity  and 
density  fluctuations.  In  Figure  28,  the  velocity  fluctuations  are 
Plotted  aa  per  Equation  (61),  and  thla  plot  should  be  compared  with  that 
Of  Figure  25.  There  1,  obviously  little  difference  between  these  two 
Illustration.,  which  can  be  attributed  to  both  the  unavoidable  data 
•carter  and  the  fact  that  the  density  gradient,  extant  in  this  measure¬ 
ment  are  too  small  to  Illuminate  the  difference.  In  any  event,  the  form 
Of  the  function  f(„)  1.  quit,  clear  and  is,  ..  expected,  very  close 
numerically  to  Its  form  at  low  sp.eds  (see  Reference  1,  Figure  7  4 
Reference  3.  Figure  9,  and  Reference  5,  Figure  17),  Including  the  Nation 
and  value  of  the  maximum..  This  form,  a  strong  argument  that  the  DEH 
correctly  describes  the  compressible  wake. 

The  effect  of  proper  correlation  1,  the  more  evident  In  the  density 
fluctuation,  of  Figure  29.  Her,  the  absence  of  large  total  temperature 
fluctuation,  collapse,  the  data  to  a  form  very  much  the  same  a,  1„ 

Figure  28,  exactly  a.  Equations  (61)  and  (71)  predict;  more  Importantly,  ' 
the  peak  In  the  den.lty  fluctuation,  shifts  In  magnitude  and  location  so 
•a  to  euke  that  posalble.  To  Illustrate  this  shift,  Figure  30  has  beer. 

_  '  * 

‘Comparison  with  low-speed  data  1.  handicapped  by  the  nonuniform  wanner 
by  which  the  wake  radius  ha.  been  normalised  In  the  referent.  ,  quoted. 


FIGURE  29.  APPARENT  SELF-PRESERVING  FORM  OF  THE  RMS  DENSITY  FLUCTUATIONS  NORMALIZED 
ACCORDING  TO  THE  DYNAMIC  EQUILIBRIUM  HYPOTHESIS. 


FIGURE  30.  EFFECT  OF  FROFKR  COORDINATE  TRANSFORMATION  ON  THE  RADIAL  FLUCTUATION  DISTRIBUTION 
AT  A  TYPICAL  AXIAL  POSITION.  IN  THE  TRANSFORMATION  SUGGESTED  BY  THE  DEH 
(BOTTOM)  THE  DISTRIBUTION  MAXIMA  IN  THE  VELOCITY  AND  DENSITY  COME  CLOSER  IN 
MAGNITUDE  AND  POSITION.  ;/»■ 


plotted.  Here  the  velocity  and  density  fluctuations  have  been  again 
plotted  as  in  Figures  25,  26,  28,  and  29;  the  shift  in  the  peak  is 
clear.  In  Figure  30,  as  in  the  previous  illustrations,  it  is  also  clear 
that  at  higher  T}  (near  the  wake  edge),  the  distribution  f  differs 
between  the  velocity  and  density  fluctuations.  This  is  attributable 
to  the  effect  of  the  intermittency . 


According  to  Equation  (67),  the  axis  rms  velocity  fluctuations, 
normalized  by  the  free-stream  velocity  u,*,  should  decay  as  the  (-2/3) 
power  of  X  only  if  multiplied  by  the  density  factor  (1-r)1^2;  such  a 
plot  is  shown  in  Figure  31.  The  agreement  with  the  predicted  decay  law 
is  unsatisfactory;  with  the  possible  exception  of  the  last  few  experi¬ 
mental  points,  the  data  decay  faster  than  the  (-2/3)  power  of  X, 
although,  in  fairness  to  the  DEH,  not  as  fast  as  the Au ^0^/Uoo  plot  of 
Figure  23  indicates. 

Better  agreement  obtains  when  the  form  of  Equation  (59)  is  tested  as 
shown  in  Figure  32.  In  spite  of  the  scatter,  the  velocity  fluctuations 
normalized  approach  clearly  the  limlt'^jg(o)  ■  0.38  which  the  DEH, 
fortified  by  Townsend's  data  (Reference  1),  indicate.  On  the  same  figure, 
the  density  fluctuations  attain  the  same  numerical  limit  as  predicted 
by  Equation  (79).  Thij  behavior  is  quite  gratifying  because  the  ordinates 
here  are  formed  by  grouping  various  independent  measured  experimental 
quantities . 


A  further  manipulation  of  the  data  of  Figure  32  removes  the  need  for  a 
numerical  choice  for  g(o)  and  reduces  the  scatter:  such  is  the  ratio  of 
the  axis  density  to  the  velocity  fluctuations  in  each  case  normalized 
with  their  mean  differences.  This  ratio  should,  according  to  Equation 
(80),  Increase  in  the  similar  region  and  tend  toward  l/S*1.18  at  very 
large  X.  The  solid  line  in  Figure  33  is  the  quantity  (l-r)(l-w)/S  as 
obtained  from  the  mean  measurements,  and  the  agreement  with  the  data  is 
very  good. 


FIGURE  32. 
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The  decay  o£  Che  wake  temperature  fluctuation  analogous  to  the  velocity 
fluctuation  decay  of  Figure  22  has  been  predicted  on  the  basis  of  the 
strong  Reynoldsjinalogy  (Equation  (50))  and  the  relevant  ratio  Is 
plotted  versus  X  In  Figure  34.  Conal.tent  with  the  small  total  temper- 
ature  fluctuations  found,  this  ratio  tends  toward  unit  (i.e.,  toward 
thejterlflcatlou  of  the  analogy  and  of  Equation  (50))  at  large  values 
of  X.  Note  that  In  this  and  several  previous  Illustrations,  a  relaxation 
distance  of  about  40  wake  (drag)  diameters  Is  found,  consistent  with 
the  mean-flow  results  (Reference  12).  other  forms  of  the  temperature 
fluctuation  have  not  been  plotted  because  the  conversion  of,  say, 
AT(o)/'T(o)  to  AT(o)/Tw  is  a  simple  matter. 

In  order  to  test  the  analogy  in  the  radial  direction,  the  ratio  of 
Figure  34  is  also  plotted  versus  T)  in  Figure  35.  There  are,  in  fact, 
two  general  types  of  behavior  this  ratio  follows:  nearer  the  base  it 
peaks  slightly  towards  unity  about  one-half  the  wake  radius  away  from 
the  axis;  farther  downstream,  the  ratio  i.  very  close  to  unity  throughout 
except  near  the  wake  edge,  where  it  increases  greatly  beyond  unity.  Both 
these  pehnomena  can  be  attributed  to  the  intermittent  nature  of  the 
turbulent  front.  In  the  first  instance  (Figure  35- a),  the  local  total 
temperature  fluctuations  ATq  are  important  enough  to  destroy  the  analogy 
within  the  turbulence  itself,  and  the  ratio  lies  sensibly  below  unity. 

Near  the  actual  wake  edge,*  however,  that  pseudoturbulent  component  due 
to  the  intermittency  Itself  emphasizes  the  density  fluctuations  with  the 
result  that  a  peak  appears  in  the  ratio.  Farther  downstream  (Figure  35-b) 
the  influence  of  ATq  has  decreased,  but  pseudoturbulent  density  fluctua¬ 
tions  due  to  the  intermittency  appear  at  the  wake  edge,  thickening  the 

arent  turbulence  diameter  (Figures  28  and  29)  and  the  ratio  is  thereby 
diiven  to  high  value. 


♦This  is  the  edge  as  defined  properly  by  intermittency  measurements, 
which  will  be  described  In  a  later  report. 
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SECTION  6 

INTEGRAL  SCALES 


The  integral  scales  of  axial  velocity  and  temperature  (density)  were 
computed  by  the  WEB-VII  program  (Appendix  E)  as  a  function  of  position, 
according  to  the  scheme  of  Equations  (47)  and  (48).  The  results  are 
shown  in  Figures  36  through  39.  As  already  discussed  in  Paragraph  4.1, 
there  are  intuitive  arguments  that  the  integral  scales  (i.e.,  the 
autocorrelation  macroscales)  should  be  locally  proportional  to  the 
characteristic  length,  L,  of  the  wake  and  thus  grow  with  the  1/3  power  of 
axial  distance.  To  date,  however,  no  information  has  existed  on  the 
numerical  differences  between  the  density  (temperature)  macroscale, A$, 
and  the  axial  velocity  macroscale,  A. ror  thelr  radl«l  variation. 

Figure  36  shows  the  variation  of  the  axis  values  of  the  A $  normalized 
alternately  with^/c^A  and  with  L.  There  is  little  doubt  that  on  the  axis 
the  ratio  A_/L  is  constant  in  the  similar  region  of  the  wake, 
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FIGURE  36,  AXIAL  VARIATION  OF  LONGITUDINAL  DENSITY  FLUCTUATION  INTEGRAL  AHTOCORRFt_att™  c 
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FIGURE  37. 
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RADIAL  VARIATION  OF  THE  LONGITUDINAL  DENSITY  FLUCTUATION 
INTEGRAL  AUTOCORRELATION  SCALE,  NORMALIZED  WITH  THE  AXIS 
VALUE  IN  THE  RELAXING  (TOF)  AND  THE  SELF-PRESERVING  WAKE  (BOTTOM) 
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RADIAL  VARIATION  OF  THE  LONGITUDINAL  VELOCITY  FLUCTUATION 
INTEGRAL  AUTOCORRELATION  SCALE,  NORMALIZED  WITH  THE  AXIS  VALUE 
IN  THE  RELAXING  (TOP)  AND  THE  SELF -PRESERVING  WAKE  (BOTTOM), 
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Th*  radial  variation  of  A.  1.  shovn  1„  Flgur.  37.  Immedlat.ly  dovnssrsam 
Of  th.  transition  tons,  the  density  nu.crosc.le  decreases  to  a  minimum  off 
tie  axis  and.  In  fact,  attains  a  value  a.  small  .s  0.5  A (T)  -  o) .  Beyond 
that  the  scale  Increase,  systematically  with  distance,  and  In  the  similar 
portion  (X  >  55.9)  It  Increases  linearly  from  the  axl.  to  the  apparent 

edge  of  the  wake  vhere  Tig  ■  2A  s  <T)  -  o).  In  term,  of  the  transverse 

scale,  L  this  implies 

AS  ^  *  °)  *  0.5  L 

As  01  -  3.5)  -  L 

with  a  linear  variation  in  between. 

As  seen  in  Figure  38,  the  axis  variation  of  the  axial  velocity  auto¬ 
correlation  macroscale At  cannot  be  pictured  effectively  because  of 
excessive  data  scatter  due  in  Urge  part  to  scatter  in  the  computer 
output  in  the  velocity  spectra  at  very  low  frequencies;  the  trend  with 
increasing  X  is  fairly  obvious,  however,  with  the  limiting  value  of  AT 
lying  between  0.5  L  and  L.  Ihere  is  less  doubt  on  the  radial  variation 

of  AT,  and  from  Figure  39  we  see  that  here  in  fact  its  value  is  practically 
constant  at  A ^  ■  AT(o)  across  the  wake. 

The  conclusion  drawn  is  that  the  expected  intimate  connection  between 
As  (and  At)  with  L  is  indeed  verified.  It  is  also  evident  from  Figure  36, 
and  to  some  extent  from  Figure  38)  that  the  equilibrium  value,  of  the 
scales  are  arrived  at  by  decrease  from  the  higher  values  near  the  transition 
region.  Ihia  is  exactly  what  is  expected  in  events  where  large-scale 
Phenomena,  such  as  laminar  instabilities,  break  up  into  the  smaller  scale 
turbulent  structure.  Very  much  che  same  results  have  been  lately  obtained 
by  Gaviglio  and  Favre  (Reference  19). 


Figure  40  shows  the  ratio  of  the  macroscale  Ag  to  the  macrbscale  Ar  In 
accordance  with  Figure.  37  and  39,  this  varies  from  about  0.5  on  the  axis 
to  about  unity  at  the  edge  of  the  wake. 


An  Important  byproduct  of  the  scale  measurement  Is  the  reassurance  that  the 
experimental  conditions  are  capable  of  producing  </alid  statistical  averages. 
The  validity  condition  obviously  Is  that  the  hot-wire  samples  a  sufficiently 
large  number  of  macroscale-size  eddies  during  Its  exposure  to  the  local 
turbulent  flow.  In  this  experiment,  such  local  conditions  are  (for  exped¬ 
iency)  measured  by  traversing  the  wire  radially  across  the  wake  (cf.  Figure  7) 
and  then  by  processing  the  data  once  all  the  necessary  traverses  are  taken. 

The  question  therefore  arises  as  to  whether  the  traverse  speed  is  so  fast 
that  the  wire  does  not  remain  at  any  particular  point  in  the  wake  to  sample 
a  sufficient  number  of  eddies.  Each  traverse  lasts  approximately  10  seconds, 
implying  that  the  wire  remains  approximately  1  second  in  a  radial  segment 
of  the  size  of  the  macroscale  (about  1  mm).  For  the  prevailing  axial  mean 
velocity  of,  typically,  50,000  centimeters  per  second,  there  are  thus  about 
a  half-million  eddies  in  the  aggregate  from  which  the  average  is  drawn  -- 
a  good  statistical  average. 


SECTION  7 

SPECTRAL-MODAL  ANALYSIS 


We  shall  now  present  the  results  of  the  spectral -modal  analysis  wherein  the 
distribution  of  fluctuation  intensities  will  be  further  expanded  into  its 
Fourier  contributions.  Before  discussing  these  results,  we  will  utilize 
the  typical  datum  from  the  VEB-VII  program  output,  shown  in  Figure  41,  in 
order  to  define  terms,  dimensions  and  units. 


There  are  three  main  representations  of  interest  for  each  mode  and  at  each 
point  in  the  wake.  First,  we  want  to  know  how  much  of  the  total  hot-wire 
signal  is  due  to  density  and  how  much  is  due  to  axial  velocity  fluctuations. 
Thus,  within  a  1-kilocycle-per-second  passband  around  each  frequency,  we 
define  byo(f)  tho  root -mean -square  density  fluctuation  contribution 
(normalized  with  the  local  mean  density)  and  byr(f)  the  corresponding 
contribution  of  the  velocity  fluctuation.  Referring  to  Figure  41,  these 
results  are  shown  on  the  second  and  third  columns  as  a  function  of  the 
frequency,  shown  in  the  first  column.  The  units  of  these  quantities  are 
(keps)  1;  to  convert  these  to  (cps)"1,  i.e.,  seconds,  one  divides  by  1000, 
and  this  has  been  done  in  the  illustrations  plotting  these  quantities.  The 


\ 


\ 


)  « 


1 


?  ”  2  *“  2f  5?5*  5?  I?  "*  or*  *  °  **“  •  0  **  *>  <•  o®  *>  «  «  <v  A.  £>  m  A. *:  <i>  o 

_  ,*  ^nrf(io«N«Ko**  »-f 

'  ••»)>»>  fvrv^.  —  r*  +  r>  r\ 


0 

•ft 

0 

O 

0 

O 

♦ 

4 

0 

9*. 

c 

4 

rv 

IT 

JS 

w- 

w- 

•4. 

S' 

« 

« 

O 

0 

< 

«-  •— 

*  *-4. 

w 

r-< 

a 

►"  H 

CM  O 

W  (3 

C.  4C- 

1  O 

O  O 

u*  • 

•  ♦ 

.0  0 

O  J\ 

n 

m  h 

•  0 

10  0 

<  0 

H>  Hi 

0 

H>  O 

O 

O  HI 

• 

•  • 

c 

0  r> 

£ 

CO 

c 

<v  t- 

*— 4 

3 

► 

4  H 

HI 

►*  Ci 

►- 

a> 

y. 

.'  at 

3  < 

► 

Cft 

0 

H  O 

0 

O  O 

0 

O  O 

♦ 

•  ♦ 

1* 

90  tO 

fO 

fv  0 

0 

CM  tf\ 

0 

f4i  H> 

•0 

0  0 

*> 

*  m 

• 

•  • 

0 

0  0 

CO 

0  -* 

H  O 

CM 

W  «• 

H  *4 

*/»  ui 

X 

*-♦  x 

13 

»“  M 

H  < 

«•  or 

»  O 

h  l 

9 

40  O 

X  X 

Z 

3  Ht 

U)  or 

in  _• 

V 

4 

■*  O 

H  O 

O 

O  O 

O 

O  O 

♦ 

•  ♦ 

O 

CM  .V 

O 

©  O 

O 

O  fO 

O 

CM  ¥> 

O 

O  CM 

Jb 

N  -O 

• 

•  • 

O 

O  O 

CM 

fi 

in  ;n 

%0 

Ul 

<  H 

X 

0  a 

H 

ft  a 

<T 

x  r 

a 

-1  j 

<m 

H  O 

0 

O  O 

0 

O  0 

CO 

•  ♦ 

• 

H  O 

* 

cm  r* 

0  ** 

#- 

0  rv 

tf>  H 

a 

O  M 

F— 

•  • 

OOOOOOOOOO  OOOOO  OOOOOOOOOO  OOOOO 
M°?T°T°CC?COOOEC  oooooo  cococ  C  000 

■‘£“i5cooooc’c'oel',®0,«,*,*0~*«>oo=ooao 

n  e  n.  ■*•  ti 

«(VOOHNH)«l(\|fl<|*H»'*li)N%4o* 

'7  *  *  O'  X> 


‘O  O  O  O  o  O  OOOOO  OOOO  o 


hhhh«4HmHHccoouo  00  o  c-  o  c  no  00  „•  c  oc  n 

O  O  O  O  O  O  O  o  *  q  q  ^  ^  ^  ^  ^  | ^ 

CJUC000000003000  OO  O  O  C  C  OOOO'JCO"  O 

*“  •  •  •  1  *  •  ••  •  ♦♦♦♦♦♦♦♦♦♦♦»♦♦♦«,  ►♦♦♦ 

•  5«  «  «  «>  <  -  ■•  ^ 


tni^«BOIVnil\xHNAINNNNI4l4n«  tr«tNt>G  MI 

ICOOOrlHHriOHnxi  «  tfl  «  ^  o  90m  C  oncim  <->  ir\  -t  —t 
I  CV  *■»  «  «  «  ^  O  O  ..  ~  ^  W  ,4  ,*  ^  ^  X  W  «*  M  £  K)  n  *  *  S  £  •  £ 


OC.OOOOOCJOO  OOOOO  OOOOOO 


OOOOO  0000 


«4v4«Hv4v*HI*4*4f4v4  *-40000  OOOOOO  OOOO  «4  H 

^oocooooooccooeoooooonoooc  cc 

\000000000,3000000  ooooooooooo  » 


C  (  C.  O  C  C  C 


v  o  c  tj  c  r. 


>v4Hff4HplHnHHC>OOOCOCC  O  U  r  O  H  r  r(  H 

>00000000000  OOOOOOOOOO  o  000 


Xoocccoooooooooooocoocnnroir  c*  *  r  r> 

JJ*  ♦♦♦♦•♦♦♦♦♦•♦♦♦♦♦•♦♦♦♦♦»♦♦  r  ,  , 
WNCVJMMmnvoiCiaAilXIlSeKiKfVnNC  l>*  t>  4>  «l  -)  rs  o«  « 


..  _  -  -  - -  —  »  f>  N  f  (V  n  O  <M»)  N  o  t  « 

^  0  •  C  O  h  h  f(  /ma  ixx  hn  4,  oorvo»i>  hoa  «or<>  r«  ** 
0r,£v<  «ftOCC*\nHCOO»  «  40  (\i  <f  fv  «  «  K  ri«  K  c  O  K 


•  OOOOOOOO  OO  00004 


>000  D  5  O  3  O 


i  ?  «  «  «  «  S!  SS  2  S!  ~  ~  ^  ^  ^  ~  w  ^  ~  ^  ^  ~  ~  ~  ^  ^ ♦ 

W IOOOOOOOOOOOOOOO oooooooooooooo 

•  O  OOOOOOOOOO  OOOOO  oooooaaann  _  —  — .  — 


V  <  OOOOO' 


•*♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 7 

3l9»|(VOCNM«(V9  otinMif)^  f  Ht  ONr  «t  ^P»  t  * 

^  *•  H  O  ♦  ©  Oi  OHO©©  HOH©©©  o  •< 

^  ?  ??  ®  ^  ^  ^  ^  ^  W  (M  O  HS  OO  fs  ©OK  ^  O ©  ©  tO 


WTF7Pjmo  HN  o«n  ©ON  «  O©  IT 

•  oooonooooon oooooo  oo 


030000 30 000003000000 
OOOOOOOO  000000000000 
OOOOOOOO  OOOOOO  OO  ooaa 


#  u  h* 

®  M 
Ul  ►*  <4 

*  *  o 


*  5  HO©WO©WC^©  wj  CW  O  H  W©?o!t  cSoS 

*  *  W  HMN^t  f^O  00S  SScM^yO  *  ~5«*#fo  ®o 

E  ffe''  •"*  HHHHfiHHa(«(VfrrHt  © 


* 


^r>JfcK 


•v 


..  :  1.  -•  a.’;.-- 


SS-V; 

a£3i 


FIGURE  41.  TYPICAL  OUTPUT  FROM  THE  VEB-VII  PR< 


representation  of  fundamental  importance  gives  the  fraction  of  the 
OU1  fluctuation  magnitude  which  1.  contained  in  the  passband  at  hand. 

“  1  ’ '  ratl°  °f  th'  -“-•»»  density  fluctuation  magnl- 

-  (  )  n  the  1-kllocycle -per -second  passband  around  the  frequency  f 

°  C°“1  ^r*1uency- Integra  ted)  me.n-.qu.re  density  fluctuation  °2  at 

t  e  same  point  in  the  wake;  note  that  the  connection  between  off)  ,„d  c 
(and  also  T(f)  andr)  have  been  given  by  Equations  (42)  through  (46) 

This  ratio,  and  the  corresponding  ratio  T2(f)/T2,  .re  f„„„d  ln  the  slxt„ 

and  seventh  columns  of  Figure  41.  Again,  the  units  here  are  In  (kcp,)'1 
or  milliseconds. 

One  can  employ  the  third  representation,  shown  in  the  eighth  and  ninth 
columns,  to  compare  the  spectra  by  the  expedient  of  normalising  the  spectral 
densities  by  their  respective  sero-frequency  limits.  These  normalised 

IZrXluTT'’  dl"en,l0nl*’S  a"d  «•  footed  as  2(f)/o2(0) 

(  )/T  (0).  In  this  form  the  spectra  can  be  compared  with  the  custom¬ 
er,  one-d Intensions  1  spectral  density  description,  such  as  shown  In 
Reference  18.  To  this  end,  the  last  two  column,  of  Figure  41  also  shiv 

appropriate  nondlmenslonallzatlon.  of  the  frequency  Into  the  form, 
f  Am 


u 
f  A 


(90) 

ns  *  <91> 

where  AT  and  As  are  the  respective  Integral  scales  and  „  is  the  local  mean 
velocity. 

With  this  introduction,  the  numerous  Illustration,  shown  of  spectral  densi¬ 
ties  are  self-explanatory.  I„  drawing  comparison,  and  conclusion,  it  should 
be  notea  that,  ln  the  axial  direction,  we  should  normally  look  for  an  evolu¬ 
tionary  adjustment  of  the  spectra  within  the  relaxation  sone  (X  <  40)  and 
perhaps  Invariant  behavior  In  the  rnglon  of  sUt-preservationT  In  the 
radial  directly,  we  draw  guidelines  from  the  results  of  the  modal  analysis 
and  the  unpublished  results  of  the  Intermittent,  measurements.  Thu.  we  are 
Interested  ln  the  following  regions;  (1)  the  wake  axis  ffl-0),  where 


f 


( 


I 


Ui 


I 
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fluid  is  fully  turbulent,  (2)  the  region  of  maximum  shear  (7)  -  l.l)  which 
coincides  with,  and  obviously  gives  rise  to,  the  maxima  in  the  fluctuation 
magnitudes  (compare  Figures  28  and  29) ,  (3)  the  location,  in  the  statistical 
sense,  of  the  wake  front  (r)  -  2.11),  and  (4)  the  wake  edge  which  can  be 
conveniently  defined  as  the  point  at  which  the  mean  velocity  distribution 
u  is  1  percent  of  its  axis  value  (i.e.,  0.01)  and  which  is  thus  located 
(Reference  12)  at  T]  -  3.28. 


A  good  first  look  of  the  distribution  of  spectral  densities  is  given  by 
Figures  42  and  43.  These  serve  to  bring  out  the  unusual  feature  of  a 
fluctuation  concentration  (peak)  at  a  scale  which,  from  Figure  42,  is  about 
6  times  and,  from  Figure  43,  about  12  times  the  integral  scales.  In 
Figure  42,  drawn  from  the  relaxation  zone  (X  *  26.5)  and  in  Figure  43, 
taken  from  the  self-preserving  region  (X  -  56.0),  this  peak  is  generally 
absent  from  the  axis  and  the  wake  edge  and  seems  to  become  most  intense  in 
the  region  1  <  T)  <  2,  i.e.,  between  the  maximum  shear  and  the  front  location. 
We  shall  return  to  the  discussion  of  this  feature  in  the  comments  which 
follow. 


Turning  to  the  spectra  on  the  axis.  Figure  44  shows  that,  initially,  the 
wire  signal  is  predominated  by  velocity  fluctuations  uniformly  greater  than 
the  density  fluctuations  at  all  frequenc  farther  downstream,  of  course, 
the  situation  is  reversed  (as  the  modal  a^»  ysis  has  already  shown), 
although  at  X  -  61.9  the  preponderance  ot  the  density  fluctuations  is  mainly 
due  to  the  larger  eddies.  In  the  same  vein,  consider  the  results  at 
X  -  26.5  of  Figure  45.  Here  OaT,  but  the  velocity  predominates  at  the 
lower  frequencies,  with  the  density  taking  over  a  range  of  higher  frequencies. 
Figure  46  shows  that  the  density  fluctuation  density  (in  terms  of  fraction 
of  the  total  local  density  fluctuation)  is  distributed  with  frequency  in  a 
way  independent  of  axial  position,  and  this  similarity  becomes  quite 
remarkable,  in  fact,  when  the  spectra  are  normalized  in  the  manner  of 
Figure  47.  By  contrau  to  the  latter,  which  presents  data  beyond  X  -  56, 
consider  that  the  density  spectral  shape  is  much  more  scattered  in  the 


FIGURE  44.  AXIAL  EVOLUTION  07  THE  LONGITUDINAL  SPECTRA  OF  RMS  DENSITY  AND  VELOCITY 
FLUCTUATIONS  (NORMALIZED  WITH  THE  LOCAL  MEAN  VALUES)  ON  THE  WAKE  AXIS. 
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LONGITUDINAL  SrZCTJU  OF  TH1  KEAN -SQUARE  DCNSITT  FLUCTUATIONS  (NORMALIZED  fc’XTB 
THE  INTEGRATED  FLUCTIA  'IONS)  ON  THE  AXIS  OP  THE  SELF -PEE  SERVING  WAKE. 
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relaxing  wake,  as  the  data  in  14.7  <  X  <  50  plotted  in  Figure  48  show. 
Unfortunately,  no  such  conclusions  can  be  drawn  from  Figure  49  for  the 
for  the  velocity  fluctuations,  due  to  the  large  scatter.* 

In  the  region  of  maximum  shear,  Figure  50  again  reflects  the  predominance 
of  density  fluctuations,  whether  early  (X  =  26.5)  or  late  (X  *  85.4)  in  the 
wake.  At  this  radial  position  in  the  wake,  the  observed  peak  sets  in  at 
full  strength.  In  Figure  51,  the  fractional  velocity  fluctuation  density 
clearly  shows  the  evolution  of  the  peak  at  X  =  26.5,  50,  and  85.4. 

Initially,  the  peak  is  located  at  higher  frequency  (n^  se  0.2)  and  occupies 
a  higher  percentage  of  the  total  fluctuation;  later,  it  shifts  to  lower 
frequencies  and  decreases  in  magnitude.  A  similar  result  is  seen  in 
Figure  52  for  the  density  fluctuations  as  they  develop  in  the  relaxation 
zone.  By  contrast,  similar  results  farther  along  in  the  wake,  shown  in 
Figure  53,  show  little  difference  from  X  to  X,  and  in'  fact  the  percentage 
fluctuation  occupied  by  the  peak  remains  quite  constant.  On  the  other 
hand,  the  nondimens ional  plot  of  Figure  54  destroys  this  constancy;  even 
so,  the  peak  is  well  discernible.  Figure  53  should  be  contrasted  with 
Figure  55  which  once  more  points  out  the  unsettled  nature  of  the  spectrum 
during  the  relaxation  process.  Inspection  of  Figure  55  reveals  the  evolution 
of  the  peak  itself  and  also  the  accompanying  variation  of  the  high-frequency 
components  (also  shown  in  Figure  52).  It  is  evident  that  Fourier  contri¬ 
butions  of  this  peak  are  to  be  found  within  a  large  range  of  frequencies:  as 
the  peak  shifts  to  the  left,  so  do  the  higher-frequency  components.  As  we 
turn  to  Figures  56  and  57,  we  observe  that  not  only  has  the  scatter  in  the 
velocity  spectra  decreased  (compare  Figure  49) ,  but  the  peak  behaves 
differently.  Specifically,  it  seems  to  disappear  completely  at  the  far 


*As  is  easily  seen  from  the  spectra,  the  velocity  fluctuation  density  results 
are  always  accompanied  by  large  scatter,  partly  because  these  fluctuations, 
being  smaller  than  those  of  the  density,  infringe  on  the  error  margin  of 
the  measurements. 
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FIGURE  52.  AXIAL  EVOLUTION  OF  THE  LONGITUDINAL  SPECTRUM  OF  THE  ME AN -SQUARE  DENSITY 
FLUCTUATIONS  (NORMALIZED  WITH  THE  INTEGRATED  FLUCTUATIONS)  IN  THE  REGION 
OF  MAXIMUM  SHEAR. 
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FIGURE  33.  THE  LONGITUDINAL  SPECTRA  OP  MEAN-SQUARE  DENSITY  FLUCTUATIONS  (NORMALIZED 
WITH  THE  INTEGRATED  FLUCTUATION)  IN  THE  REGION  OP  MAXIMUM  SHEAR  OF  THE 
SELF -PRESERVING  WARE. 
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FIGURE  55.  THE  LONGITUDINAL  SPECTRA  OF  THE  MEAN-SQUARE  DENSITY  FLUCTUATIONS  (NORMALIZED 
AT  ZERO  FREQUENCY)  IN  THE  MAXDIM-SHEAR  REGION  OF  THE  RELAXING  WAKE.  NOTE 
DEVELOPMENT  OF  PEAK. 
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FIGURE  57.  THE  LONGITUDINAL  SPECTRA  OP  THE  MEAN-SQUARE  VELOCITY  FLUCTUATIONS 
(NORMALIZED  AT  ZERO  FREQUENCY)  IN  THE  MAXIMUM-SHEAR  REGION  OP  THE 
SELT -PRESERVING  WAKE. 
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distances  (X  -  67.7,  85.4,  91.3).  We 
mechanise  is  responsible  for  this  peak 
density  rather  than  the  velocity. 


therefore  conclude  that,  whatever 
,  it  operates  mainly  through  the 


A  generally  identical  sort  of  behavior  of  the  spectral  densities  can  be 
found  along  the  turbulent  front  (TJ  =  2.1),  as  can  be  seen  from  the  following 
group  of  illustrations.  Of  these,  Figure  58  shows  the  well-established 
decreases  of  the  densities  in  terms  of  fractions  of  the  local  mean  values, 
and  the  tendency  of  this  decrease  to  occur  primarily  in  the  larger  eddies! 
Figures  59,  60,  and  61  use  the  alternate  spectral  representations  to 
illuminate  the  evolution  of  the  spectra  with  axial  distance.  Again,  the 
shift  of  the  peak  to  lower  frequencies  and  its  decrease  with  distance  are 
perfectly  clear.  It  is  furthermore  obvious  that,  as  noticed  previously, 
the  density  fluctuation  density  retains  the  evidence  of  this  peak  for  a 
perceptibly  longer  axial  distance. 


As  a  final  region  of  interest,  consider  the  wake  edge  which,  in  view  of  our 
previous  findings,  can  be  thought  as  lying  around  T\  =  3.3.  The  nondimensional 
spectra  of  the  density  fluctuations  for  this  radial  position  have  been 
plotted  in  Figure  62  in  the  self-preserving  region.  This  plot  verifies 
what  was  already  indicated,  to  some  extent,  by  Figures  42  and  43,  i.e., 

that  the  spectra  have  returned  to  a  normal  shape,  the  peak  being  hardly 
discernible . 


To  formalize  what  is  learned  about  the  radial  distribution  of  spectral 
densities,  we  present  Figure  63,  where  the  nondimensional  density  spectra 
have  been  plotted  for  all  radial  positions,  at  the  farthest  downstream 

position  investigated  (at  X  -  91.3).  The  appearance  of  the  peak  here  is 
quite  weak. 


An  overall  conclusion  which  can  be  drawn  from  the  inspection  of  the  spectra 
is  that,  the  observed  peak  aside,  there  is  little  difference  between  them 
and  what  is  usually  observed  in  low-speed  turbulent  flows  in  which  a  certain 
amount  of  homogeneity  and  isotropy  are  obtained.  In  Figure  64,  we  compare 
the  present  results  (typically  those  shown  in  Figure  63)  with  spectra 
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FICUKI  63.  LONGITUDINAL  SPECTRA  OF  THE  MEAN-SQUARE  DENSITY  FLUCTUATIONS  (NORMALIZED 
WITH  THE  INTEGRATED  FLUCTUATIONS)  IN  THE  SELF-FRESERVINC  WAKE. 


obtained  by  Favre  (Reference  20)  behind  a  grid.  A  commonly  ueed  analytical 
expression  for  the  one-dimensional  velocity  fluctuation  spectrum  is  also 
shown  for  comparison  These  three  spectra  differ  little  for  eddy  sizes 
larger  than  the  integral  scale,  but  there  is  considerable  disagreement 
among  the  three  at  the  higher  frequencies.  In  fact,  the  fluctuation 
density  thus  nomalized  shows  that  in  the  compressible  wake  studiec  the 
fluctuation  density  is,  for  eddies  the  size  of  the  integral  scale,  one-half 
as  large  us  at  low  speeds. 

The  new  and,  to  a  certain  degree  unexpected  phenomenon,  is  t lie  appearance 
of  the  peak  or  prominence  in  the  spectra.  We  have  seen  that  tills  phenomenon 
shows  certain  organized  traits  which  will  be  used,  in  the  ensuing  discussion, 
to  discover  its  origin. 

We  have  noted  that  the  peak  behaves  as  follows:  (1)  it  appears  strongest 
about  halfway  between  the  wake  axis  and  the  wake  edge,  diminishing  toward 
the  axis  and  the  edge;  (2)  it  appears  most  markedly  in  the  initial 
(relaxation  stage)  and  thereafter  decays;  (3)  the  peak  frequency  is  initially 
quite  high  (tending  toward  a  length  scale  approaching  the  Integral  scale 
of  the  turbulence)  but  thereafter  decreases;  (4)  the  peak  is  more  pronounced, 
although  not  by  very  much,  in  the  density  rather  than  the  velocity  spectrum. 
To  these  we  add  two  additional  observations:  (3)  when  the  peak  frequency 
is  plotted  versus  the  radius  at  various  X  as  in  Figure  65,  we  note  that  in 
the  usual  dimensional  units  the  frequency  is  on  the  order  of  50  kilocycles 
per  second  at  the  beginning,  and  that  its  previously  noted  decrease  with  X 
brings  it  closer  to  an  average  value  of  between  20  and  30  kilocycles  per 
second.  From  Figure  65  one  further  sees  that  the  value  of  the  peak 
frequency  tends  to  attain  a  maximum  at  about  the  midpoint  between  the 
axis  (T)  ■  o)  and  the  edge.  (We  shall  return  to  the  us*  of  these  phenomena.) 
Last,  the  hot-wire  signal  was  passed  through  an  lntermittency  circuit  whose 
function  war  to  rectify  the  turbulence  component  and  mark,  by  an  on-off 
signal,  the  entrance  and  exit  of  the  wire  into  the  zones  of  turbulence. 

The  spectrum  of  this  output  signal  was  measured  and  was  found  to  be 


FIGURE  65.  PEAK  SPECTRUM  FREQUENCY  VARIATION  WITH  WAKE  4.0IUS. 
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characteristically  identical  with  the  total  wire  output.  The  peak, 
especially,  survived  the  rectification  process  and  appeared  intact  in  the 
rectified  spectrum.  The  latter  is,  in  fact,  the  same  as  the  Schmidt-trigger 
output  spectrum  whic1'  Corrsin  and  Kistler  (Reference  21)  measured  in  the 
turbulent  interface  of  a  jet  and  found  to  be  of  the  Poisson  type. 

The  latter  test  provides  the  strongest  first  clue  as  to  the  origin  of  the 
peak,  namely  that  it  is  due  to'the  fluctuations  arising  from  the  corrugations 
of  the  interface  rather  than  the  turbulence  within  these  corrugations.  If, 
roughly  speaking,  the  intermittent  signal  is  thought  oi  as  a  superposition 
of  turbulence  on  a  random  train  of  rectangular  waves,  then  the  superposition 
will  carry  over  to  the  spectrum  which  should,  in  principle,  be  similarly 
decomposable  into  the  bona-fide  spectrum  and  the  pseudoturbulent  front 
signal.  Returning  to  point  (5)  and  Figure  65,  we  add  here  the  observation, 
that  when  the  "zeros"  of  the  rectified  wave  were  measured,  they  yielded  a 
characteristic  frequency  (crossings  of  the  wire  into  the  turbulent  front) 
also  of  about  30  kilocycles  per  second,  in  agreement  with  the  spectrum 
measurement  results.  This  frequency  naturally  also  decreased  near  the 
axis  and  again  near  the  wake  edge,  implying  that  Figure  65  shows  the 
behavior  of  the  front  rather  than  of  the  turbulence  itself. 

Comments  (1)  through  (4)  can  now  be  put  into  proper  context.  If  our 
hypothesis  is  correct,  then  the  front  pseudoturbulence  should  be  most 
intense  in  regions  of  the  most  intense  intermittency,  which  occurs  at  about 
the  middle  of  the  wake  radius;  this  is  in  agreement  with  Comment  (1).  The 
pseudoturbulence  decay,  commented  on  in  (2)  is,  of  course,  fully  predictable. 
Observation  (4)  is  also  easy  to  interpret  because  the  momentum  exchange 
between  the  large  corrugations  of  the  front  and  the  external  flow  is 
probably  much  more  efficient  than  the  exchange  of  heat.  On  the  other  hand,  - 
the  exact  interpretation  of  the  frequency  shift  is  in  some  doubt.  If  our 
hypothesis  is  true  that  the  front  (interface)  is  indeed  weakly  periodic, 
then  this  frequency  should  be  characteristic  of  the  wavelength  of  the  front 
corrugations.  Because  the  axial  velocity  of  the  front,  itself  unknown, 


-118- 


must  ar>  j  a  the  axial  direction,  at  Least  initially,  there  is  little  useful 
information  the  initial  peak  frequency  shift  can  supply.  At  the  same  time, 
the  apparent  relaxation  of  the  peak  frequency  location  to  about  30  kilocycles 
per  second  (see  Figure  65)  is  consistent  with  observations  (Reference  22) 
of  the  frequency  of  "zero  occurrences"  which  was  also  of  about  30  kilocycles 
per  second.  As  we  saw  from  the  spectra,  this  corresponds  to  a  scale  about 
12  times  higher  than  the  typical  macroscale  size,  which  implies  a  Length  of 
between  6L  and  12L  where  L  is  the  transverse  scale  of  the  wake.  From 
measurements  of  the  intermittency  and  zeros,  the  former  agrees  with  the 
microscale  of  the  wake  front  and  is  thus  advanced  here  as  a  representative 
value  of  the  wavelength  of  the  front. 

In  summarizing,  the  only  new  phenomenon  emerging  from  the  spectral  measure¬ 
ments  is  not  necessarily  the  inevitable  mixing  of  the  turbulent  signals  with 
the  pseudoturbulence  of  the  front,  but  rather  that  the  latter  is  not, 
apparently,  truly  random.  This  organization  of  the  Latter  into  a  weakly 
periodic  structure  serves  at  least  as  an  identification  of  the  pseudo- 
turbulent  action  and  could  be  of  potential  use  in  separating  the  turbulent 
from  the  pseudoturbulent  phenomena. 

The  discussion  of  the  modal -spectral  analysis  will  be  concluded  with  a 
brief  glance  of  the  spectral  distribution  of  the  velocity-temperature 
cross-correlation  coefficient,  as  presented  in  Figure  66.  It  is  seen  that 
anti-correlation  is  generally  most  perfect  for  the  largest  eddies  lying 
on  the  wake  axis  near  the  transition  zone,  while  it  is  poorest  for  the 
smallest  eddies,  again  on  the  axis,  but  far  downstream.  It  is  especially 
interesting  that  correlation  is  generally  best  for  the  largest  eddies  at 
any  rate,  verifying  Morkovin's  conclusion  that  these  eddies  are  the  ones 
responsible  for  transport  across  the  gradient. 

The  errors  involved  in  this  measurement  are  in  general  very  difficult  to 
estimate  accurately.  Estimates  can  only  be  made  of  the  sources  of  error, 
and  the  specific  techniques  reducing  such  errors  can  be  equally  well 
identified.  Thus,  major  error-producing  mechanisms  were  (1)  the  inherent 
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FIGURE  66.  SPECTRAL  RESOLUTION  OF  THE  VELOCITY -TEMPERATURE  CROSS -CORRELATION  COEFFICIENT. 


alignment  difficulties  of  axi- symmetric  models,  (2)  slight  differences  in 
re-installing  and  re-aligning  the  mode’  each  time  it  had  to  be  removed  from 
the  test  section,  (3)  not  fully  satisfactory  wind  tunnel  flow  (from  the 
standpoint  of  uniformity),  (4)  differences  in  wire  characteristics,  since 
the  inevitable  wire  breakages  required  the  uses  of  many  wires,*  and 
(5)  unavoidable  imperfections  in  the  transfer-function  technique.  On  the 
other  hand,  the  care  exercised  in  calibrating  wires  and  identifying  the 

strain  gage  problem  has  been  of  great  help  in  reducing  anomalies  in  the 
data. 

By  neglecting  the  sound  fluctuations  in  analyzing  the  data  we  have  in 
effect  stated  that  the  normalized  pressure  fluctuations  are  insignificant 
compared  to  the  density  and  temperature  fluctuations.  In  this  context  the 
self-consistency  of  the  experiment  can  be  checked  by  computing  the  pressure 
fluctuations  induced  by  velocity  fluctuations  as 


f  *  -  1  »2  (*)’ 


where  M  is  the  Mach  number.  The  following  tabulation  presents  results 
obtained  for  two  radial  locations  at  X  -  14.7.  The  agreement  with  the 
imposed  condition  seems  very  satisfactory,  particularly  away  from  the  axis. 


Comparison  of  Pressure  and  Density  Fluctuations 


(ftu/u) 

0.24 

0.055 


(Ap/p) 

0.075 

0.016 


MLE1 

0.120 


The  large  fluctuation  levels  encountered  naturally  cast  some  doubt  on  the 
validity  of  the  technique  and  the  analysis  used  to  reduce  the  data.  Speci¬ 
fically,  at  least  three  possible  sources  of  error  arise:  (1)  the  wire 
heat  transfer  characteristic  may  be  highly  nonlinear  over  the  range  of  the 

*  ^  ULthat  Wlth  crltlcal  or  delicate  model  geometries 

deiJc*teLhot  wlr«*»  obtains  the  data  as  fast  as  practical 
without  disturbing  the  model  In  the  meantime.  P 
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fluctuations,  (2)  the  wire  may  become  sensitive  to  the  transverse  as  well 
as  the  axial  component  of  the  fluctuations  in  the  velocity  vector,  and 
(J)  the  fluctuation  modes  interact  nonlinearly  among  themselves  so  that 
the  modal  analysis  into  entropy  and  vorticity  fluctuations  is  invalid. 
Because  higher-order  moments  in  the  turbulence  variables  greatly  complicate 
the  formulas,  an  exact  computation  of  the  errors  due  to  these  nonlinear¬ 
ities  is  impossible,  although  some  rough  estimates  can  be  made.  Thus,  the 
heat  transfer  characteristic  of  the  wire  is  more  linear  (Re0  ~  Nu0)  at  the 
lower  Reynolds  numbers  than  at  the  higher  ones  and  thus  derivatives  such 
as  d2NuG/dRe2  are  very  small;  the  instrument  is  therefore  more  linear  at 
these  low  Reynolds  numbers.  The  inclusion  of  nonaxial  velocity  fluctua¬ 
tions,  roughly  speaking,  sets  an  error  in  (Au/u)  of  about  15  percent  for 
fluctuation  on  the  order  of  0.2,  but  the  error  decreases  rapidly  with  the 
fluctuation  level,  i.e.  ,  as  one  moves  away  (axially  or  radially)  from  the 
point  of  maximum  level.  The  nonlinear  interaction  of  modes  is  difficult 
even  to  estimate  except  perhaps  for  interaction  involving  sound  (see 
Reference  15). 

There  is,  too,  the  problem  of  possible  distortion  of  the  mean  measurements 
because  of  the  high  turbulence  level.  Fortunately,  the  reduction  of  the 
mean  data  rested  on  the  static  pressure  which  is  largely  insensitive  to 
the  fluctuations  and  on  the  total  temperature  which  was  not  found  to 
fluctuate  much.  These  facts,  coupled  with  the  insensitivity  of  the  pitot 
tube  to  the  turbulence,  indicate  that  the  (p^.,  p,  TQ)  group  of  measurements 
should  be  always  preferred  over  the  (pT,  p,  Rt0)  group  for  mapping  the 
mean  field  in  highly  turbulent  flows. 
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SECTION  8 


CONCLUSIONS 


On  the  basis  of  the  work  performed  to  date,  the  following  conclusions  can 
be  drawn  regarding  the  turbulent  wake  investigated: 

(1)  The  transition  to  turbulence  of  the  axi- 
symmetric  was  clearly  triggered  by  the 
amplification  of  a  laminar  instability. 

The  highest  local  fluctuation  levels 
observed  lay  in  or  immediately  downstream 
of  the  transition  zone. 

(2)  Th^  overwhelming  majority  of  the  turbulence 
properties  tended  to  relax  into  a  self¬ 
preserving  form  within  a  region  of  about 

40  virtual  diameters  downstream  of  transition. 

However,  the  possibility  exists  that  tho 
relaxation  distance  is  in  general  a  function 
of  the  free-stream  flow  parameters. 


L4 


-  \ 

v  '  *  k  ;  . 
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(3)  It  has  been  verified  that,  on  account  of  the 
small  total  temperature  fluctuations  measured, 
the  so-called  strong  Reynolds  analogy  holds 
for  wakes  produced  by  adiabatic  bodies. 

The  radial  distributions  of  the  temperature 
(density)  and  axial  velocity  fluctuations, 
when  scaled  by  factors  drawn  from  the  Reynolds 
analogy  and  the  dynamic  equilibrium  hypothesis, 
appear  to  be  identical  with  each  other  and  to 
their  low-speed  counterparts.  This  similarity, 
occurring  in  the  physical  rather  than  the 
Howarth-Dorodnitzyn  plane,  verifies  that  the 
equilibration  nrocess  is  local  and  maintained 
by  the  larger  eddies.  Thus,  the  velocity 
and  temperature  fluctuations  are  well  anti¬ 
correlated,  mainly  for  the  large  eddies. 

The  axis  values  of  tht  velocity  and  temperature 
fluctuations  are  likewise  in  agreement  with 
the  qualitative  and  quantitative  features  of 
the  dynamic  equilibrium  hypothesis. 

(6)  The  velocity  auto-correlation  macroscale  i a, 
within  the  scatter,  equal  to  the  transverse 
wake  scale,  L;  it  is  fairly  constant  with 
radius.  In  the  self-preserving  region  (by 
contrast),  the  density  auto-correlation 
macroscale  is  about  0.5  L  on  the  axis  and 
increases  towards  L  at  the  wake  edge;  these 
two  scales  are  not  therefore  everywhere 
equal  to  each  other.  By  way  of  comparison, 
the  maximua  shear  zone  occurs  at  1.1  L,  the 
statistical  wake  front  position  is  at  2.1  L, 


(5) 


(7)  The  normalized  spectral  densities  of  the 
velocity  and  temperature  are  virtually 
Identical  and  uniform  with  radius  in  the 
self-preserving  region.  They  are  identical 
with  the  low-speed  wake  velocity  fluctuation 
spectrum  at  the  low  frequencies,  but  of  de¬ 
creased  intensity  at  frequencies  corresponding 
to  (and  higher  than)  the  macroscale. 

(8)  The  spectral  densities  art  distorted  by  the 
pseudoturbulence  generated  by  the  wake  front 
intermittency ,  which  is  weakly  periodic  with 
a  wavelength  of  about  6  L.  The  relaxation 
length  of  this  pseudoturbulence  appears  to 

be  longer  than  that  for  the  turbulence  within 
this  front. 

It  should  be  reemphasized  that  the  study  presented  here  deals  with  an 
essentially  steady-state  system  where  Eulerian  statistics  are  performed  by 
an  instrument  stationary  in  a  flowing  gas.  By  this  process,  the  resulting 
map  of  statistical  properties  contains  a  mixture  of  both  the  turbulence 
and  the  pseudoturbuleuce  of  the  front.  This  composite  picture  is  both 
legitimate  and  convenient  for  practical  applications  but  it  invites  further 
work  by  which  the  turbulence  due  to  the  front  will  be  separated  from  the 
bona-fide  turbulence  inside  it.  In  this  new,  equally  relevant  representation, 
where  the  intermittency  factor  figures  prominently,  the  radial  distributions 
of  properties  reported  herein  will  be  altered  by  subtracting  the  front 
pseudoturbulence  and  adding  an  amount  of  turbulence  normally  unaccounted 
because  of  intermittency.  Such  work  wiT.l  be  reported  at  a  later  date. 
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APPENDIX  A 

THE  TRANSFER  FUNCTIONS 


A.  1  HOT  WIRE 

The  transfer  function  W  of  the  hot-wire  la  given  by 


e(f) 


wire 


e(f) 


flow 


W(f) 


-  [  1  +  A  TT2  f2  M2 


-1/2 


;  M  s  time  constant. 


(A.  1) 


where  the  phase  shift  is  not  included;  here  f  is  the  frequency  and  M  the 
wire  time  constant.  The  latter  Is  given  by 
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Here  C  Is  the  overall  heat-capacity  of  the  wire 
C  -  (density)  x  (specific  heat)  x  (volume) 

°rRr  the  wen-calibration  resistance-temperature  slope,  Rtf  the  heated  and 
Rr  the  equilibrium  wire  resistance,  I  the  current,  A^  the  overheat 
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and  e  the  finite  circuit  factor  which  was  measured  to  be  0.008  in  our  case. 
The  quantity  af  is  zero  in  our  case  and  thus  we  can  re-write  (A. 2)  as 


2 

cpTir  l  A' 

M  -  - 1  seconds  (A. 4) 

r  r  I 

Where  the  specific  heat  c,  density  p,  wire  radius  r,  and  length  l  etc.,  are 
given  in  cgs  units.  The  transfer  function  W  is  then  given  by  combining 
(A. 1)  with  (A. 4). 


At  this  Juncture,  it  is  appropriate  also  to  give  the  limiting  value  of  M  at 
zero  current.  By  (A. 3)  and  (A.*)  we  get 
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(A. 5) 


where  e  is  the  dc  voltage.  We  recall  that  at  zero  current  the  Nusselt 
number  is 
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and  since  at  I-*0 


and 
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the  combination  (A. 5)  and  (A. 6)  gives 


A.  2  AMPLIFIER 

The  zero-frequency  gain  G0  of  the  amplifier  is  given  in  Figure  A-l.  The 
amplifier  transfer  function  G/G0  has  been  calibrated  for  different  attenua¬ 
tion  settings,  is  shown  on  Figure  A-2,  and  has  been  curve-fitted  by 

-  [l  +  d-25)  10”12  f2+  (0.75)  10‘17  f3]’1  (A. 8) 

where  f  is  in  cycles  per  second. 

A. 3  COMPENSATOR 

When  the  compensator  time-constant  dial  is  set  at  Ms  (which  is  equal  to 
Mi  of  program  WEB-IT.)  the  compensator  transfer  function  is 
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1  +  (0.705)  10"2  f  Ms  +  (0.369)  10'3  f2  M2  -  (3.2)  lo'6  f3 
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which  has  been  measured  to  500  kilo-ycles  per  second  and  is  shown  in 
Figure  A-3. 
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(A.  9) 


A. 4  RMS  METER 

The  transfer  function  of  the  Ballantine  rms  meter  set  at  300  millivolt  full 
scale  is 
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(A. 10) 


where  E  is  the  dc  voltage  output  from  the  mean  square  Ja:k  on  the  meter  in 
millivolts.  For  the  1-volt  full-scale  setting  this  is 

_  _  _ 

*VTVM  0A3  0.685  J7T94  (A.  11) 
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These  curves  are  shown  on  Figure  A-4. 
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FIGURE  A-2.  TRANSFER  FUNCTION  OF  THE  AMPLIFIER, 
OF  ZERO - FREQUENCY  GAIN, 


SHOWN  TO  BE  INDEPENDENT 
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A. 5  WAVE  ANALYZER 
A. 5.1  BANDWIDTH 

As  Figure  A-5  shows,  checks  made  on  the  frequency-independent,  selectable 
bandwidth  of  the  wave  analyzer  snowed  excellent  filtering  characteristics 
at  0.2-,  3-  and  1-kilocycle-per -second  width;  the  latter  is  used  for  the 
measurements . 

A . 5 . 2  TRANSFER  FUNCT ION 

The  transfer  function  H,  shown  in  Figure  A-6,  was  found  to  be  independent 
of  frequency;  tests  arr  always  made  at  the  settings  ABSOLUTE,  NORMAL,  1  VOLT 
f.s.,  -10  decibels  down  and  using  a  1000-ohm  load.  In  this  case,  an  rms 
input  of  e  millivolts  is  converted  to  dc  millivolts  V  by 

V 

e  ’  H 
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APPENDIX  B 


THE  WEB-II  COMPUTER  PROGRAM 


The  general  purpose  of  this  program  is  to  nullify  the  frequency-response 
deficiencies  of  the  apparatus  in  measuring  the  frequency- integrated  voltage 
output  of  the  hot-wire  at  each  current  and  each  point  in  the  waks.  This  is 
done  as  follows  in  each  case:  the  output  spectrum  is  recorded,  the  overall 
transfer  function  measured,  and  an  error-ratio  formed  from  the  spectrum 
r<easurement  alone,  indicating  how  much  the  measured  integrated  wire  output 
should  be  Increased  to  account  for  apparatus  deficiencies. 

Each  spectrum,  taken  in  the  range  of  0  to  500  kilocycles  per  second  is 

first  broken  into  7-kilocycle-|>er-sf  cond  intervals  (the  latter  Increase  at 

the  higher  frequency  end)  ana  the  analyzer  signal  and  noise  output  noted 

and  the  spectral  density  of  the  signal  determined;  integrated  over  the 

2 

spectrum  this  gives  the  integrated  wire  output  In  the  meantime, 

the  transfer  functions  of  the  amplifier  and  compensator  are  evaluated  at 
each  frequency.  The  transfer  function  of  the  wire  is  found  as  follows: 
the  wire  time  constant  outside  the  wake  is  known  (for  each  current  and 
axial  position),  as  is  its  current -voltage  curve.  The  local  time  constant 
at  each  point  in  the  wake  is  then  found  by  the  formula  given  it.  the  main 

taxt,  and  the  local  transfer  function  of  the  wire  is  computed  and  combined 

WB  ■'■<1  \  |D  ■  i-.'L  jk 

with  the  others  to  form  the  overall  transfer  function  CTT .  The  value  of 


I 


the  latter  at  each  frequency  is  now  combined  with  the  net  wave  analyzer 
output  at  that  frequency  tc  give  the  ldP2:  hot-wire  output  at  that  fre¬ 
quency;  Integration  then  leads  to  the  Ideal  integrated  wire  output  e*2 

_  VTVM 

at  that  point  and  current.  The  error  ratio  J  Is  formed  by  e  and 
2  o  VTVM 

**VTVM  and  combined  wlth  the  wlre  output  eQyTyM  measured  with  the  rms 

voltmeter  to  correct  the  latter. 

At  each  current  and  point  In  the  wake  the  program  output  presents,  among 

other  quantities,  the  wire  time  constant  M  (M(MSEC)),  the  ideal  and  actual 

integrated  wire  outputs  e*JTVM  (VTVM*2)  and  ej^  (VTVM2),  respectively, 

as  computed  from  the  spectrum,  the  same  quantity  as  measured  directly 
2 

eOVTVM  *  the  error  ratio  J,  and  the  final  corrected  integrated  wire 

output  e2  (E2).  Further,  at  each  station  and  current  the  following  are 

among  quantities  given  as  a  function  of  the  frequency  f  (FKCPS):  the  net 
2 

analyzer  output  V  (f),  the  wire,  amplifier,  and  compensator  transfer 

—  functions  W,  G/Gq  and  C,  respectively,  the  overall  transfer  function  OTF, 

2 

and  the  net  wire  signal  power  density  in  (volts)  /cps. 

As  noted  from  the  above  remarks,  this  program  provides  a  meticulous  and 
detailed  view  of  the  behavior  of  the  complete  electromechanical  system 
over  the  entire  matrix  of  operating  situations.  In  addition  to  the  new 
insight  it  offers  into  the  technique  of  anemometry,  it  affords  numerous 
cross  checks  of  the  procedure  and  increases  accuracy.  More  importantly, 
the  output  can  now  serve  as  a  starting  input  for  a  combined  modal-spectral 
analysis  (i.e.,  the  spectra  of  the  velocity  or  density  fluctuations  alone), 
a  task  normally  so  cumbersome  that  it  has  been  attempted  only  once  previously. 

The  procedural  diagram  of  WEB-II  is  shown  on  Figure  B-l. 


APPENDIX  C 

THE  WEB-III  COMPUTER  PROGRAM 


The  purpose  of  this  program  wss  to  process  the  turbulence  data  through  the 
response-restoration  process  of  the  WEB-II  program  and  to  prepare  them  for 
the  eventual  analysis  Into  frequency-integrated  modes  (performed  ».<lth  the 
WEB- IV  program). 


The  basic  Inputs  to  this  program  are  the  15  radial  traces  of  integrated 
raw  hot-wire  turbulence  signal  (1  for  each  of  15  currents)  and  the  corres¬ 
ponding  15  traces  of  mean  wire  voltage;  the  error  ratio  J  is  also  put  In. 

The  coordinates  of  each  trace  (distinguished  by  T  for  the  turbulence  and 
S  for  the  mean  or  steady  traces)  are  supplied  by  an  analog-to-digital  card- 
punch  as  "counts;"  by  appropriate  conversion  factors  (F^  and  Fg),  these 
traces  are  reduced  into  curves  of  voltage  ET  versus  radial  distance  Y. 
Utilizing  the  transfer  function  of  the  rms  voltmeter,  the  root-mean-square 
of  the  turbulent  signal  C31ST)  is  obtained.  This  process  is  performed  by 

a  taking  discrete  points,  20  in  all,  for  each  x  station.  After  subtraction 
of  the  noise  RMSN,  the  actual  integrated  wire  voltage  output  at  these 
points  is  known,  and  can  now  be  corrected  for  the  ideal  response  by  using 
the  error  ratio  J  (cf.  WEB-II  prograsi).  As  one  recalls,  J  had  been  obtained 
at  8  currents  I,  not  15  as  available  in  this  program,  and  also  at  10,  rather 


than  20  radial  positions;  interpolation  thus  became  necessary,  and  the  J's 
so  obtained  were  then  used  to  get  the  ideal  integrated  turbulence  signal. 


I 


In  similar  fashion,  the  mean  wire  voltage  was  obtained  for  each  point 
using  level  L  from  the  plotter  curves,  corrected  into  e  to  account  for  the 
line  resistance  (1.04  ohm),  squared  and  combined  with  e^  to  form  the  ideal 
nondimens ional  frequency-integrated  (but  amplified)  wire  turbulence  output 
at  each  current  and  point;  this  was  further  reduced  by  dividing  by  the 

zero-frequency  amplifier  gain  G  .  Finally,  the  wire  overheat  parameter  A', 

2.-22°  w 
against  which  (e  /e  )/G^  is  given  in  the  output  was  formed. 

The  computing  procedure  for  this  program  is  shown  in  Figure  C-l 
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FIGURE  C 


1.  THE  WEB-III  PROGRAM 


APPENDIX  D 

THE  WEB-IV  PROGRAM 


The  purpose  of  this  program  was  to  supply  the  modal  analysis  of  the  radial 
distribution  of  frequency-integrated  turbulent  fluctuations;  final  objec¬ 
tive  was  therefore  to  give  the  velocity  and  density  fluctuations  and  their 
correlation  as  a  function  of  radius  and  axial  position  along  the  wake. 


G 


As  a  first  step,  the  program  prepared  a  list  of  the  flow  properties  at  each 
point  in  the  wake,  as  derived  by  the  WEB-I  program;  the  location  of  these 
points  (20  per  x  station)  were  the  same  as  in  WEB-III.  Various  combina¬ 
tions  of  these  flow  properties  were  computed  preparatory  to  insertion  in 
the  wire  sensitivity  coefficients.  Since  the  Reynolds  number  ReQ  and  the 
Re  -Nuq  anu  ReQ  -  T)  relations  were  known  fr  i  the  flow  calibration,  the 
logarithmic  derivatives  (Reo/T))  (ffldReo)  and  (Reo/Nuo)  ONuo/dReo)  were 
computed,  as  was  the  wire  temperature  Ty  at  each  current  and  the  attendant 
values  of  T  and  K.  The  circuit  parameter  Cf  was  computed  from  the  pre- 
viously  measured  A\  At  this  point,  sufficient  information  existed  for 
computing  the  wire  sensitivities  to  the  mass  flux  (em)  and  total  tempera¬ 
ture  (eT)  fluctuations.  Use  of  the  local  Mach  number  M(Y)  also  allowed 
computation  of  aisiilar  aenaitivitieo  e0  and  eT  to  the  entropy  (temperature) 
and  vorticity  (velocity*  '(fluctuations,  respectively. 
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2  2-2 

At  this  point,  the  net  ideal  integrated  mean  square  wire  voltage  e  /GQe 
from  WEB-III  was  introduced  at  each  point  and  current,  giving  15  second- 
degree  polynomials  (for  the  mass  flux  and  representation)  with  known 
coefficients,  and  a  similar  group  for  the  entropy-vorticity  representation. 
A  least-squares  solution  o c  these  did  not  work  out,  and  therefore  they  were 
solved  by  fitting  them  to  a  second-degree  curve,  obtaining  the  mass  flux 
and  total  temperature  fluctuations  and  their  correlation  m,  t,  and  r^t , 
respectively,  the  corresponding  vorticity  (T)  and  entropy  (a)  fluctuations, 
and  the  correlation  rT(J.  Finally,  these  fluctuations  were  also  given 
normalized  to  the  local  velocity  and  density  differences  (u«  -  u(0))  and 
-  p(0)  across  the  wake,  as  put  in  from  the  WEB-I  program. 

The  computations  procedure  for  this  program  is  shown  in  Figure  D-l. 
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APPENDIX  E 


THE  WEB-VII  COMPUTER  PROGRAM 


Tho.  purpose  of  the  WEB-VII  computer  program  Is  (1)  to  perform  the  combined 
modal-spectral  analysis  of  the  turbulent  flow  and  (2)  to  fur thernnanipu late 
the  modal  results  of  the  WEB-IV  program  in  order  to  get  parameters  interest 
ing  from  a  theoretical  viewpoint.  Thus,  this  program  supplies  spectra, 
macroscales,  turbulent  energy  ratios  and  the  like. 


E.l  MODAL- SPECTRAL  ANALYSIS 

Here  we  want  to  find  the  fluctuation  density  at  each  frequency  and  at  each 
point  in  the  wake.  The  main  input  is  the  corrected  "iW  voltage  fluctua¬ 
tion  G2e2,  a  function  of  the  frequency,  which  is  given  in  millivolts 
sqcarec  per  kilocycle  and  appears  on  the  last  column  of  the  WEB- II  outputs. 
This  quantity,  it  should  be  notad,  appears  multiplied  by  the  amplifier 
zero -frequency  gain  (squared)  G2,  and  is  first  divided  by  the  latter  to^ 
give  EZFLOW  (mv)2/kc.  By  dividing  the  latter  by  (IMA)*  (ROHMNET)  (mv) 
we  get  the  normalized  (modally  integrated)  fluctuation  demity  (per  kc) 
at  the  particular  frequency  Z2,  and  by  further  dividing  by  the  entrcpy 
fluctuation  coefficient  squared  (ES)2  (obtained  by  interpolation  for  the 
particular  current  IMA  and  YSTATION)  we  get  the  Kovasznay-Morkovin  variable 
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Now  we  can  proceed  with  the  modal  analysis  at  each  frequency  by  solving 
simultaneously  the  eight  equations 

Z2  -  D62  +  E6  +  F 

ni  ni  ni 

where  i-1,  2 . 8  corresponds  to  each  of  the  8  wire  currents  employed  in 

WEB-II  for  that  frequency;  here,  of  course,  the  variable  6n^  is  the  ratio 
of  the  sensitivity  coefficients 


This  solution  gives  the  results 


*  .  E 

rCT  "  27*  0* 

with  the  star  denoting  spectral  density  -  e.g.,  T*  is  the  rms  velocity 

fluctuation  (normalized  to  the  mean  local  velocity  at  that  point)  per 

2  ^2 

kilocycle.  Corollary  quantities  formed  by  this  program  are  T*  and  <7*  and 


The  latter  are  the  mean-square  velocity  and  density  (temperature)  fluctua¬ 
tion  density,  normalized  with  the  square  of  the  local  mean  velocity  and 
density,  respectively,  per  cycle  per  second.  One  quantity  of  interest, 
put  out  by  WEB- VII,  is  tha  comparison  of  tha  quantities  T*/o+.  Another 
is  the  spectral  density  of  each  mode  squar-sd  compared  with  the  frequency- 
integrated  mode 


t*2  „  a*2 

- -  end  - - 
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where 


where  the  1000-factor  appears  In  order  to  convert  to  integrand  into  cycles 
per  second  (since  T*2  and  0*2  ire  obtained  within  kilocycle  "windows"). 

The  latter  two  integrals  should  obviously  be  equal  to  the  integrated  vel¬ 
ocity  and  density  fluctuations,  respectively,  at  that  point  in  the  wake. 

The  pertinent  integral  scales  (macroscales)  are  computed  from  the  relations 
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where  u  is  the  local  flow  velocity  (obtained  from  the  WEB-I  output)  and 
s2(0)  and  t2(0)  r.re  the  limit  values  of  s2  and  t2  at  zero  frequency.  In 
the  present  instance  these  quantities  were  evaluated  at  7  kilocycles  per 
second.  Because  both  numerator  and  denominator  in  the  latter  two  equa¬ 
tions  are  each  normalized  to  the  local  mean  values,  the  definition  of  the 
scales  is  independent  of  these  values. 

2  2 

This  program  finally  outputed  the  nondimen sional  spectra  s  /s  (0)  and 
0)  versus  the  nondimens ional  frequencies 
P(CPS)  At 


F(CPS)  Ag 


u 


Note  that  the  latter  outputs  present  best  comparison  between  the  spectra 
of  these  two  modes,  whereas  the  output  represents  best  the  portions  of  the 


spectrum  contributing  most  heavily  to  each  mode. 


Other  outputs  gave  comparisons  of  the  scales  with  the  transverse  wake 

1/2 

scale  L  (from  WEB-V)  and  the  virtual  wake  (drag)  diameter  (C^A) 

This  program  also  provided  for  the  manipulation  of  certain  quantities 
put  out  by  the  WEB-IV  program;  the  following  were  formed  as  a  function  of 
position  in  the  wake: 


(1) 


The  quantity  TAU/SIGMA  and  the  related  quantity 


SRA 


SIGMA  1 

TAU  CY-1)M2 


where  M  is  the  local  Mach  number. 


(2)  The  rms  velocity  and  density  fluctuations  at 
each  XSTATION,  each  normalized  with  the  axis 
value.  (See  UPRIME/UPRIMO  and  RHOPRIM/RHOPRIMO. ) 


(3) 


The  local  average  turbulent  energy  divided  by  the 
local  mean  kinetic  energy 


~  s  ENERGY  - 
Em 


3T2  -  3.46  raTTO 


(4) 


where  we  have  assumed  Isotropy  in  order  to 
compute  the  turbulent  velocity  vector  from  the 
measured  longitudinal  component,  and  where 
third-  and  higher-order  moments  are  neglected. 

Other  quantities  manipulating  the  fluctuations 
according  to  the  dynamic  equilibrium  hypothesis. 


Typical  outputs  appear  in  Figure  41. 


